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PREFACE 


The University of Texas at San Antonio (UTSA) is a relatively new university. It 
was established in 1969 and opened for classes in 1973. As the only comprehensive public 
university serving the nation's ninth largest city, it was and is vital to San Antonio and the 
entire South Texas Region. In 1982, just eight years ago, an undergraduate engineering 
program was established at UTSA with the support of the community and its leaders 
Today, all three undergraduate engineering programs are ABETaceredited and serve about 
800 students, a significant percentage of whom are Hispanic. The future includes a new 
engineering budding, containing laboratory facilities and equipment, planned to open m 
January, 1991. Furthermore, a graduate program has just been put m place at the M.S. 
level and one is planned at the Ph.D. level. The first Master's Degree students enrolled in 

Fall, 1989. 

Naturally, the engineering research environment is just developing at UTSA. Now, 
thanks in great measure to the UT System support and this ongoing NAS A grant, good 
progress is being made. Specifically, the purchase of a UT System Cray-X-MP in March, 
1986 and a second one in December, 1988 has provided a world-class analytical and 
numerical research environment not ordinarily available to a new university. As a result the 
UTSA Supercomputer Network Research Facility (SNRF) was developed by the pnncipa 
investigator, Dr. Lola Boyce. This has aUowed the successful completion of this research 
project, an early one of its kind at UTSA. 

This NASA research grant has allowed three Mechanical Engineering students, 
Thomas Lovelace, Callie (Scheldt) Bast and Eddie Aponte, to work directly with the 
principal investigator. Dr. Boyce, providing them with a quality research experience they 
would otherwise probably not have had. All students have expressed an interest in 
continuing their education at the graduate level. 

In conclusion, and in view of the significant accomplishments in fundamental 
research, enhancement of the engineering research environment at UTSA, and direct 
support of Mechanical Engineering students, it is hoped that the proposed extension of t is 
grant will receive favorable consideration at NASA. The principal investigator sincere y 
thanks NASA for funding this second year grant. 
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ABSTRACT 


This report presents the results of & second year effort of a research program 
conducted for NASA-LeRC by The University of Texas at San Antonio (UTSA). The 
research included development of methodology that provides probabilistic lifetime strength 
of aerospace materials via computational simulation. A probabilistic phenomenological 
constitutive relationship, in the form of a randomized multifactor interaction equation, is 
postulated for strength degradation of structural components of aerospace propulsion 
systems subjected to a number of effects or primitive variables. These primitive variables 
often originate in the environment and may include stress from loading, temperature, 
chemical or radiation attack. Time may also interact with them, producing cr ee p and 
fatigue. In most cases, strength is reduced as a result This multifactor interaction 
constitutive equation is included in the computer program, PROMISS. Also included in the 
research is the development of methodology to calibrate the above-described constitutive 
equation using actual experimental materials data together with a multiple linear regression 
of that data, thereby predicting values for the empirical material constants for each effect or 
primitive variable. This methodology is included in the computer program, PROMISC. 
Actual experimental materials data were obtained from the open literature for a number of 
nickel-base superalloys and a few metal-matrix materials, materials typically of interest for 
aerospace propulsion system components. 
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1.0 INTRODUCTION 


This report presents the results of a second year effort of a research program 
entitled "Development of Advanced Methodologies for Probabilistic Constitutive 
Relationships of Material Strength Models, Phase 2." TTiis research is sponsored by the 
National Aeronautics and Space Administration-Lewis Research Center (NASA-LeRC). 

The principal investigator is Dr. Lola Boyce, Associate I^fessorof Mechanical 
Eneineering The University of Texas at San Antonio (UTSA). The objective of the 
research program is the development of methodology that provides probabilistic lifetime 
strength of aerospace materials via computational simulation. 

As part of this second year effort, a probabilistic phenomenological constitutive 
relationship, in the form of a randomized multifactor interaction equation, is postulated for 
strength degradation of structural components of aerospace propulsion systems subjected to 
a number of effects or primitive variables. These primitive variables often originate m the 
environment and may include stress from loading, temperature, chemical or radiation 
attack Time may also interact with them, producing creep and fatigue. In most cases, 
strength is reduced as a result. Also included in the research is the development ot 
methodology to calibrate the multifactor interaction constitutive equation using actual 
experimental materials data together with a multiple regression of that data, thereby 
predicting values for the empirical material constants for each effect or primitive variable. 
Section 2.0 summarizes the theoretical background for the research. 

The above-described randomized multifactor interaction constitutive equation is 
included in the computer program, PROMISS. Calibration of the equation by multiple 
linear regression of the data is included in the computer program, PROMISL . 1 hese 
programs were developed using the UTSA Supercomputer 'Network Rese^ch Facility 
(SNRF) Cray X-MP, completely compatible with the NASA LeRC Cray X-MP. Thus, 
these new versions of the programs will execute on the current NASA-LeRC 
supercomputer facilities. The latest versions (Ver. 2.0) of these programs are obtainable 
from the principal investigator at the address given on the cover page of this report, l o 
verify program performance, several sample problems for each program accompany this 
report (see APPENDIX 1, Section 6.0 and APPENDIX 2, Section 6.0 and the enclosed 
floppy disks). IMSL Ver. 10 [1] subroutines are utilized. 

Sections 3.0 through 8.0 address specific tasks described in the proposal for this 
research [21. Specifically, Section 3.0 presents and discusses cases for analysis that 
resulted from a sensitivity study, utilizing the PROMISS "fixed capability. The cases 
show the effect on probabilistic lifetime strength for each of several effects or pnmitive 
variables (quasi-static stress, mechanical fatigue, thermal fatigue and creep) at both high 
and low temperature conditions. Section 4.0 describes and references the complete 
computer documentation for both PROMISS and PROMISC, contained in this report as 
User Manuals (see APPENDIX 1 and APPENDIX 2). Section 5.0 describes and 
discusses the increased capability of PROMISS, Ver. 2.0, over the original version.Ver. 

1 0 written at NASA-LeRC during the Summer of 1988. Secnon 6.0 describes and 
discusses the increased capability of PROMISC, Ver. 2.0 over the onginal version, Ver. 
1.0, written at NASA-LeRC during the Summer of 1988. Section 7.0 describes and 
discusses the materials data collected from the open literature to be used in developing da a 
for the PROMISS Resident Database. Section 8.0 discusses the precautions necessary to 


*PROMISC was written and referenced as "RAMS 11 during the Summer of 1988 and in the proposal for this 
research. 


combine PROMISS and PROMISC and the future research needed to consider such a 
combined code. 

A paper was produced documenting much of the effort of this second year research 
program. This paper is entitled "Probabilistic Constitutive Relationships for Material 
Strength Degradation Models", by L. Boyce and C C Chamis. It was presented at the 
30th Structures, Structural Dynamics and Materials Conference, Mobile, AL, April, 1989 
and is published in the Proceedings. It has also been submitted to the AIAA Journal of 
Propulsion and Power. 
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2.0 THEORETICAL BACKGROUND: 

PROBABILISTIC PHENOMENOLOGICAL CONSTITUTIVE RELATIONSHIP FOR 
LIFETIME STRENGTH OF MATERIALS 


Recently a general phenomenological constitutive relationship, for composite 
materials subjected to a number of diverse effects or primitive variables has been 
postulated to predict mechanical and thermal material properties [3, 4, 5, 6]. The resulting 
multifactor interaction constitutive equations summarize composite micromechamcs theory 
and have been used to predict material properties for a urudirecuonal fiber-reinforced 
lamina, based on the corresponding properties of the constituent materials. 

These equations have been modified to predict the mechanical property of strength 
for one constituent material due to "n" diverse effectspr primitive variables. These effects 
could include both time-independent and time-integrated primitive variables, such as 
mechanical stresses subjected to both static and impact loads, thermal stresses due to 
temperature variations and thermal shock, and other effects such as chemical reaction or 
radiation attack. They might also include other time-dependent primitive variables such as 
creep, mechanical fatigue, thermal aging, thermal fatigue, or even effect such useswnal 
attack (see APPENDIX A, Primitive Variables, Symbols, and Units). For most of these 
primitive variables, strength has been observed to decrease with an increase m the variable. 

The postulated constitutive equation accounts for the degradation of strength due to 
these primitive variables. The general form of the equation is 


c n 

f=n 

*0 i=l 


Ajf' A, 


A iF‘ A io 


( 1 ) 


where A,, Aif and Aio are the current, ultimate and reference values of a particular effect, 
a; is the value of an empirical constant for the i * effect or primitive variable, n is the 
number of product terms of primitive variables in the model, and S and Sq are the current 
and reference values of material strength. Each term has the property that if the current 
value equals the ultimate value, the current strength will be zero. Also, if the current value 
equals the reference value, the term equals one and strength is not affected by that variable. 

This deterministic constitutive model may be calibrated by an appropriately curve-fitted 
least squares multiple linear regression of experimental data [7], perhaps supplemented by 
expert opinion. Ideally, experimental data giving the relationship between effects and 
strength is obtained. For example, data for just one effect could be plotted on log-log 
paper. A good fit for the data is then obtained by a linear regression analysis. This is 
illustrated schematically in Figure 1. The postulated constitutive equation, for a single 

effect, is then obtained by noting the linear relation between log S and log ^ Ap — 


Af - A J 
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Fig. 1 Schematic of experimental data illustrating 
the effect of one primitive variable on strength. 


as follows: 



( 2 ) 


Note that the above equation (2) is for a primitive variable that lowers strength. If a 
variable raises strength, the exponent is negative. 

This general constitutive model may be used to estimate the strength of an aerospace 
propulsion system component under the influence of a number of diverse effects or 
primitive variables. The probabilistic treatment of this equation includes randomizing the 
deterministic multifactor interaction constitutive equation, performing probabilistic analysis 
by simulation and generating probability density function (p.d.f.) estimates for strength 
using a non-parametric method, maximum penalized likelihood [8,9]. Integration yields 
the cumulative distribution function (c.d.f.) from which probability statements regarding 
strength may be made. This probabilistic constitutive model predicts the random strength 
of an aerospace propulsion component due to a number of diverse random effects. 

This probabilistic constitutive model is embodied in two FORTRAN programs, 
PROMISS ('Probabilistic Material Strength Simulator) and PROMISC (Effibabilistic 
Material Strength Calibrator); see Final Technical Report, APPENDIX 2. PROMISS 
calculates the random strength of an aerospace propulsion component due to as many as 
eighteen diverse random effects. Results are presented in the form of probability density 
functions and cumulative distribution functions of normalized strength, S/So- PROMISC 
calculates the values of the empirical material constants, 

PROMISS includes a relatively simple ’’fixed" model as well as a "flexible" model. 
The fixed model postulates a probabilistic constitutive equation that considers the primitive 
variables given in Table 1. The general form of this constitutive equation is given in 
equation (1), wherein there are now n = 7 product terms, one for each effect or primitive 
variable listed above. Note that since this model has seven primitive variables, each 
containing four values of the variable, it has a total of twenty-eight variables. The flexible 
model postulates a probabilistic constitutive equation that considers up to as many as n = 18 
product terms for primitive variables. These variables may be selected to utilize the theory 
and experimental data currently available for the specific strength degradation mechanisms 
of interest. The specific effects included in the flexible model are listed in Table 2. Note 
that in order to provide for future expansion and customization of the flexible model, six 
"other" effects have been provided. 


5 



Table 1 Primitive variables available in the fixed model 


i^^intitiv^^^i^tiv^™" 
Variable Variable Type 


1 Stress due to static load 

2 Temperature 

3 Chemical reaction 

4 Stress due to impact 

5 Mechanical fatigue 

6 Thermal fatigue 

7 Creep 


Table 2 Primitive variables available in the flexible model 


A. Environmental Effects 

1. Mechanical 

a. Stress 

b. Impact 

c. Other Mechanical Effect 

2. Thermal 

a. Temperature Variation 

b. Thermal Shock 

c. Other Thermal Effect 

3. Other Environmental Effects 

a. Chemical Reaction 

b. Radiation Attack 

c. Other Environmental Effect 

B. Time-Dependent Effects 

1. Mechanical 

a. Creep 

b. Mechanical Fatigue 

c. Other Mech. Time-Dep. Effect 

2. Thermal 

a. Thermal Aging 

b. Thermal Fatigue . 

c. Other Thermal Time-Dep. Effect 

3. Other Time-Dependent Effects 

a. Corrosion 

b. Seasonal Attack 

c. Other Time-Dep. Effect 
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The considerable scatter of experimental data and the lack of an exact description of 
the underlying physical processes for the combined mechanisms of fatigue, creep, 
temperature variations, and so on, make it natural, if not necessary to consider probabilisnc 
models for a strength degradation model. Therefore, the fixed and flexible models 
corresponding to equation (1) arc "randomized", and yield die "random normalized l material 
strength due to a number of diverse random effects or primitive variables. Note that for 
the fixed model, equation (1) has the following form. 


S/Sq = f(A if, At, Aio, A 2 F, A 2 , A20»— *A7F. A 7 , A 70 ) 


(3) 


where A;, Ajf and Alo are the ultimate, current and reference values of the i* of seven 
effects or primitive variables as given in Table 1. In general, this expression can be written 

as, 

S/So = f(Xi), i = 1.-. 28, (4) 


where the Xi are the twenty eight independent variables in equation (3). Thus, the fixed 
model is "randomized" by assuming all the independent variables, Xj, 1 — 1,..., 28, to be 
random and stochastically independent. For the flexible model, equation (1) has a form 
analogous to equations (3) and (4), except that there are as many as seventy-two 
independent variables. Applying probabilistic analysis to either of these randomized 
equations yields the distribution of the dependent random variable, normalized material 
strength, S/So- 

Although a number of methods of probabilistic analysis are available, [8] simulation 
was chosen for PROMISS. Simulation utilizes a theoretical sample generated by numerical 
techniques for each of the independent random variables. One value from each sample is 
substituted into the functional relationship, equation (3), and one realization of normalized 
strength, S/So, is calculated. This calculation is repeated for each value in the set of 
samples, yielding a distribution of different values for normalized strength. 

A probability distribution function is generated from these different values of 
normalized strength, using a non-parametric method, maximum penalized likelihood. 
Maximum penalized likelihood generates the p.df. estimate using the method of maximum 
likelihood together with a penalty function to smooth it [9]. Finally, integration of the 
generated p.d.f. results in the cumulative distribution function, from which probabilities 01 
normalized strength can be directly observed. 

PROMISS includes computational algorithms for both the fixed and the flexible 
probabilistic constitutive models. As described above, PROMISS randomizes the 
following equation: 


r=n 

0 i=i 


'Nf 


-A. 


^iF' 


( 5 ) 
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where 


Af'A 1 


Af 



is the i* effect, Ai, Ajf and Ajo are random variables, aj is the i* empirical material 
constant and S/So is normalized strength. There are a maximum of eighteen possible 
effects or primitive variables that may be included in the model. For the flexible model 
option, they may be chosen by the user from those in Table 2. For the fixed model option, 
the primitive variables of Table 1 are chosen. Within each primitive variable term the 
current, ultimate and reference values and the empirical material constant may be modeled 
as either deterministic (empirical, calculated by PROMISC), normal, lognormal, or Wiebull 
random variables. Simulation is used to generate a set of realizations for normalized 
random strength, S/So, from a set of realizations for primitive variables and empirical 
material constants. Maximum penalized likelihood is used to generate an estimate for the 
p.d.f. of normalized strength, from a set of realizations of normalized strength. Integration 
of the p.d.f. yields the c.d.f. Plot files are produced to plot both the p.d.f. and the c.d.f. 
PROMISS also provides information on S/So statistics (mean, variance, standard deviation 
and coefficient of variation). A resident database, for database rather than user input of 
empirical material constants, is also provided. 

PROMISC (see Final Technical Report, APPENDIX 2) performs a multiple linear 
regression on actual experimental or simulated experimental data for as many as eighteen 
effects or primitive variables, yielding regression coefficients that are the empirical material 
constants, a;, required by PROMISS. It produces the multiple linear regression of the log 
transformation of equation (3), the PROMISS equation. When transformed it becomes 



or 


log S = log So 


X * lo 8 [ ' A » ] . 

w LAiF - AioJ 


(7) 


where 


A c -A 

iF l 

Af*Ao 


a. 

. 1 


is the i* effect, Aj, Ajf and Ajo are primitive variable data and a; is the i* empirical 
material constant, or the i* regression coefficient to be predicted by PROMISC. Also, 
log So is the log transformed reference value of strength, or the intercept regression 
coefficient to be predicted by PROMISC, and log S is the log transformed strength. 


8 


Experimental data for up to eighteen possible effects, as given in Table 2, may be included. 
The primitive variable data may be either actual experimental data or expert opinion, directly 
read from input, or simulated data where expert opinion is specified as die mean and 
standard deviation of a normal or lognormal distribution. The simulated data option for 
input data was used in the early stages of code development to verify correct performance. 
The input data, whether actual or simulated, is read in and assembled into a data matrix. 
From this data matrix, a corrected sums of squares and cross products matrix is computed. 
From this sums of squares and cross products matrix, and a least squares methodology, a 
multiple linear regression is performed to calculate estimates for the empirical material 
constant, aj, and the reference strength. So- These arc the regression coefficients. 

PROMISC includes enhancements of the multiple linear regression analysis to 
screen data from "outliers" and collinearities, determine "how well" the data fit the 
regression, quantify the importance and relative importance of each factor in the postulated 
constitutive equation, eq. (1), as weU as check assumptions inherent in the use of multiple 
linear regression. Further details are provided in the Final Technical Report, Section 6.0, 
NASA Grant No. NAG 3-867, Supp. 2, ^Probabilistic Lifetime Strength of Aerospace 
Materials via Computational Simulation." 
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3.0 CASES FOR ANALYSIS BY PROMISS 


This section presents cases for analysis that show the effect on probabilistic lifetime 
strength for each of several effects or primitive variables (quasi-static stress, mechanical 
fatigue, thermal fatigue and creep) at both high and low temperature conditions. The 
resulting sensitivity study utilized the PROMISS code "fixed" capability. A base line room 
temperature problem was established utilizing values derived from expert opinion and 
material properties typical of a cast nickel-base superalloy. The input for this problem is 
given in Table 3 and is similar to that used for the PROMISS example problem reported 
earlier [10]. For this room temperature problem, quasi-static stress was varied about 5% 
above and below the 90 ksi base line value. The results of PROMISS in the form of a 
c.d.f. is given in Figure 2. Then the base line value of 90 ksi for quasi-static stress was re- 
established, while mechanical fatigue, thermal fatigue and creep were each varied 5%, 5% 
and 90%/espectively about their base line values. The PROMISS output, again in the 
form of c.d.f. plots, for these results together with Figure 2 is given in Section 12.0, 
APPENDIX 3 as Figures A.l to A.4. A second base line problem was established that 
was identical to the first, except that temperature was set to about 57% of Tm. For this 
high temperature problem, quasi-static stress was again varied about 5% above and below 
the 90 ksi base line value. The results of this PROMISS case are given in Figure 3. Then 
the base line value of 90 ksi for quasi-static stress was reestablished, while mechanical 
fatigue, thermal fatigue and creep were each varied 5%, 5% and 90%, respectively about 
their base line values. The PROMISS output for these results together with Figure 3 is also 
given in Section 12.0, APPENDIX 3 as Figures A.5 to A.8. 

Conclusions drawn from this sensitivity study are summarized below. Raising 
temperature from room temperature to 57% Tm significantly decreases lifetime strength 
regardless of the primitive variable perturbed. Also, raising strength temperature somewhat 
decreases the spread of the three c.d.f. curves regardless of the primitive variable 
perturbed. Qualitative observations such as these are appropriate from current PROMISS 
example problems, however experimental data from actual interaction experiments must be 
used to obtain values for the empirical material constants before specific quantitative 
conclusions can be drawn for particular aerospace materials. 
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Table 3 PROMISS "fixed" model base line room temperature problem. 


Effector Variable 

Primitive Variable Symbol 


Quasi- 

SSF 

Static 

<*S 

Stress 

OSO 

P 

Temperature 

Tf 

T 

TO 

q 

Chemical 

Rp 

Reaction 

R 

RO 

r 

Impact 

SDF 

<7D 

CTDO 

s 

Mechanical 

NMF 

Fatigue 

Nm 

NmO 

t 

Thermal 

NtF 

Fatigue 

Nt 

Nto 

u 

Creep 

tCF 

tc 

tco 

V 


Units Distribution 


ksi 



Lognormal 

ksi 

Lognormal 

ksi 

Lognormal 

N/A 

Normal 

oF 

Normal 

oF 

Normal 

oF 

Normal 

N/A 

Normal 

N/A 

Normal 

N/A 

Normal 

N/A 

Normal 

N/A 

Normal 

N/A 

Normal 

N/A 

Normal 

N/A 

Normal 

N/A 

Normal 

log of cycle 

Lognormal 

log of cycle 

Lognormal 

log of cycle 

Lognormal 

N/A 

Normal 

log of cycle 

Lognormal 

log of cycle 

Lognormal 

log of cycle 

Lognormal 

N/A 

Normal 

hours 

Lognormal 

hours 

Lognormal 

hours 

Lognormal 

N/A 

Normal 


Mean Standard Deviation 

(Value) (% of Mean) 


130.0 

6.5 

5.0 

90.0 

4.5 

5.0 

-2.9 

-0.145 

5.0 

0.5 

0.015 

3.0 

2732.0 

82.0 

3.0 

68.0 

2.04 

3.0 

68.0 

2.0 

3.0 

0.5 

0.015 

3.0 

1.0 

0.003 

0.3 

0.02 

0.0006 

3.0 

0.001 

0.00003 

3.0 

0.5 

0.015 

3.0 

1.0 

0.003 

0.3 

0.1 

0.003 

3.0 

0.001 

0.00003 

3.0 

l 

O 

La 

- 0.015 

3.0 

7.0 

0.7 

10.0 

3.5 

0.35 

10.0 

1.0 

0.1 

10.0 

0.5 

0.015 

3.0 

3.0 

0.3 

10.0 

2.3 

0.23 

10.0 

1.0 

0.1 

10.0 

0.5 

0.015 

3.0 

10,000.0 

500.0 

5.0 

100.0 

3.0 

3.0 

0.083 

0.0025 

3.0 

0.5 

0.015 

3.0 





Figure 2 COMPARISON OF VARIOUS LEVELS OF UNCERTAINTY OF 
STRESS ON (ksi) PROBABLE STRENGTH FOR A 
NICKEL-BASED SUPERALLOY AT ROOM TEMPERATURE 
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Figure 3 COMPARISON OF VARIOUS LEVELS OF UNCERTAINTY 
OF STRESS (ksi) ON PROBABLE STRENGTH FOR 
A NICKEL-BASED SUPERALLOY AT 0.57 Tm 


4.0 PROMISS AND PROMISC DOCUMENTATION 


Complete documentation for both PROMISS AND PROMISC has been written, 
including both theory and user's instructions. This documentation has been appended to 
this Final Technical Report as APPENDIX 1 PROMISS USER MANUAL and 
APPENDIX 2 PROMISC USER MANUAL. These user manuals have also been issued 
under separate cover. 

The manuals provide detailed instructions for input, description of output, and 
worked example problems. Example problems are given for both "fixed" and "flexible" 
models. The theory portion of the manuals contains a description of the probabilistic 
phenomenological constitutive relations for lifetime strength of materials, [1 1, 12], the 
randomized multifactor interaction constitutive equation, relevant probabilistic and statistical 
theory, a description of all effects or primitive variables available in the code, and a 
discussion of all algorithms. Manual appendices give tables listing effect or primitive 
variable names together with their symbols, FORTRAN names and units. Other 
appendices list sample problem input and output files and IMSL subroutines used by the 
codes. The PROMISS manual also gives a sample SAS/GRAPH program for plotting 
results. Both manuals also include program notes giving information helpful to the user in 
understanding the program. 
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5.0 INCREASED CAPABILITY OF PROMISS 


The capability of the initial version of PROMISS was increased beyond the initial 
version developed at NASA Lewis during Summer, 1988. This new version is known as 
PROMISS, Ver. 2.0 and is fully documented in APPENDIX 1 of this report and under 
separate cover. PROMISS, Ver. 2.0 now includes calculation of the mean, variance, 
standard deviation and coefficient of variation of normalized strength (or lifetime strength) 
and incorporation of the Weibull as a statistical distribution to describe the uncertainty of 
any or all effects or primitive variables. Also a Resident Database was included in 
PROMISS for the empirical material constants, at Thus, rather than the user selecting 
values for these material constants, he/she may elect to have the Resident Database supply a 
value. Currently all empirical material constants in the Resident Database are set to 0.5, a 
value set by expert opinion. As further materials data become available for specific 
aerospace materials, updated values will be incorporated into the Database. 
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6.0 INCREASED CAPABILITY OF PROMISC 


The capability of the initial version of PROMISC was increased beyond the initial 
version developed at NASA Lewis during Summer, 1988. This new version is known as 
PROMISC, Ver. 2.0 and is fully documented in Section 10.0, APPENDIX 2 of this report 
and under separate cover. PROMISC, Ver. 2.0 now includes enhancements of the multiple 
linear regression analysis to determine "how well" the data fit the regression, to quantify 
the importance and relative importance of each effect or primitive variable in the postulated 
multifactor interaction equation, to screen data for "outliers" and collinearities, as well as to 
check assumptions inherent in the use of multiple linear regression. 

Multiple linear regression is a least squares technique for fitting data to a linear 
equation. It can be carried out mathematically whether it is appropriate or not Therefore, 
it can be used to test the validity of the model inherent in the multifactor interaction 
equation, by calculating "how well" the data fit the multiple linear regression through an 
Analysis of Variance (ANOVA) Table. 

Through the use of the IMSL [1] subroutine RSTAT, procedures were included in 
PROMISC that, not only assess the validity of the model inherent in the multifactor 
interaction equation, but also produce diagnostics to check for violations of the 
assumptions inherent in the use of multiple linear regression. Subroutine RSTAT supplies 
an Analysis of Variance (ANOVA) Table, which quantifies the amount of variation in the 
data that can be explained by the multiple linear regression equation. This table is therefore 
an effective method for testing the model validity inherent in die multifactor interaction 
equation. Subroutine RSTAT also supplies a Sequential Statistics Table, which quantifies 
the importance of each factor in the multifactor interaction equation. RSTAT can be 
enhanced with the capability to do a stepwise regression on the multifactor interaction 
equation. The relative importance of each term is assessed through a p-value which 
appears in the Table under the Column headed by "Prob. of Larger F". The closer that 
each p-value is to zero then the more predictive capability that associated factor has. 

Likewise, data may need to be screened for aberrant observations known as 
"outliers", through the IMSL routine RCASE. One of the most effective methods for 
outlier detection is Cook's Distance which is given in the subroutine as Cook's D. The 
problems associated with collinearities in the data can be explained through principal 
components regression. 

The subroutine RSTAT provides an Inference on Coefficient Table, which 
embellishes the results of the Sequential Statistics Table with a column headed by "Variance 
Inflation Factor". This is a diagnostic procedure for checking for the violation of 
assumptions used in multiple linear regression. The Case Analysis Table (from the IMSL 
subroutine RCASE) is useful in detecting outliers (Cook's D) and collinearities in the data. 
The residuals and studentized residuals can be used to test the assumption of normality. In 
this way, one can test the validity of all the assumptions in the multiple linear regression 
equation. A thorough discussion of two examples is given in APPENDIX 2, PROMISC 
User Manual, Section 3.0. 
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7.0 MATERIALS DATA FOR RESIDENT DATABASE 


7.1 Introduction 

In order to anchor the empirical material constants, a*, in the multifactor interaction 
equation to particular aerospace materials of interest, it is necessary to collect experimental 
fate No actual experiments were planned for this project, but instead data were collected 
from the open literature. All data obtained from the literature for this project were plotted 
and displayed in Section 13.0, Appendix 4. Data for nickel-base superalloys were far more 
abundant than that for metal matrix materials. Data for INCONEL 718 were also more 
abundant than for any other superalloy and were available for three primitive effects, 
namely high-temperature tensile* , mechanical fatigue and creep. Thus, certain INCONEL 
718 data was selected and plotted in other forms, one of which was the same as that used 
by the multifactor interaction equation in PROMISS and PROMISC (see Section.7.4). 


7.2 Literature Search 

A computerized literature search of nickel-base superalloys and metal matrix 
composites was conducted to obtain existing experimental data on various material 
properties. Useful data on high-temperature tensile, mechanical fatigue and creep 
properties were found in eight of the sources (see Section 7.5 containing the material 
bibliography page). 


7.3 Materials Data Discussion 

Data on high-temperature tensile, mechanical fatigue and creep properties were 
obtained for nickel-base superalloys, however, only creep properties were found for metal 
matrix composites. Although there is evidence to indicate that fatigue and creep processes 
interact to produce a synergistic response [1], no data on combined effects were found. In 
addition, several sources provided substantial data on INCONEL 718, a nickel-base 
superalloy. Therefore, this superalloy will be discussed separately from the other nickel- 
base superalloys. All acquired data were entered into computer input files from which 
plots could be made. Plot files were created using S AS/GRAPH software [2]. Plots from 
both nickel-base superalloy and metal matrix composite data are provided in Section 13.0, 
Appendix 4. 


7.3.1 INCONEL 718 

As stated above, high-temperature tensile, mechanical fatigue and creep data were 
obtained for the nickel-base superalloy, INCONEL 718. The data resulted from tests done 
on various hot and cold worked specimens of INCONEL 718. Tests were conducted on 
sheets of INCONEL 7 1 8 and hot rolled bars of the superalloy. Some high-temperature 
tensile and creep property data resulted from tests performed on notched specimens with a 
stress concentration factor, K t , ranging from 2.3 to greater than 20. Plots of tensile, 
fatigue and creep data for INCONEL 718 are given in Section 13.2, Appendix 4. 


* also referenced as "Temperature" in other parts of this Final Technical Report and in the User Manuals. 


17 



7.3.2 Other Nickel-Base Superalloys 


High-temperature tensile data were found for the thirteen nickel-base superalloys 
given in Table 4. 


Table 4 Nickel-Base Superalloys 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 


NldCEL-BASE SUPERAli^Y 
ASTROLOY 
HASTELLOY X 
IN-102 
IN-939 

INCONEL 600 
INCONEL 601 
INCONEL X-750 
MAR-M421 
MM-006 
NIMONIC 80A 
NIMONIC 115 
PWA 1480 SC 
RENE 41 


Data on fatigue properties were foiuid for only one superalloy, INCONEL 601. As with 
the INCONEL 718 specimens, test specimens of the other nickel-base superalloys were 
subjected to prior hot and cold work as noted on the plots. These plots are provided in 
Section 13.3, Appendix 4. 


7.3.3 Metal Matrix Composites 

Metal matrix composite data were far more scarce than nickel-base superalloy data. 
Only test data on creep properties were found. Data on high-temperature tensile or fatigue 
properties could not be found for any metal matrix composite. Data on creep behavior were 
obtained for Tungsten/Niobium and Tungsten/Niobium-1 Percent Zirconium composites. 
Two types of Tungsten fibers, 218CS and ST300, were tested as reinforcement material. 
Creep behavior plots of the metal matrix composites, as well as the Tungsten fibers, are 
provided in Section 13.4, Appendix 4. 
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7.4 Data Selected for Presentation in PROMISS and PROMISC Form 


7.4.1 High-Temperature Tensile Data 

Plots of high-temperature tensile, mechanical fatigue and creep data for the nickel- 
base superalloy, INCONEL 718, are provided in this section. Figure 4 shows the effect of 
temperature on the yield strength of the material. As expected, the yield strength of the 
material decreases as the temperature increases. 

Figure 5 is a plot of the data from Figure 4 that has been normalized with respect to 
melting temperature and yield strength at room temperature. The data along the x-axis were 
normalized by dividing each temperature value by the melting temperature, Tm, of the 
superalloy The data along the y-axis were normalized by dividing the yield strength valu 
bydiematerial’s yield strength at room temperature. Note inFigineS, that the strength 
begins to drastically decrease at a temperature that is about 50% of the melting temperature. 

Figure 6 is a plot of the log of the same data (Fig. 4) with only the x-axis 
normalized. The data along the x-axis were normalized with respect to both melong 
temperature, Tm, and a reference temperature (taken as room temperature), To. In order to 
calculate the normalized values, the difference between the melting tempera ure ligand a 
reference temperature. To, was computed This value was then divided by *e diffetence 
between the melting temperature and each temperature value, T, from Fig. 4. Finally, the 
log of each normalized value was plotted along the x-axis. Thus, the x-axis values we e 

calculated from , . 

l(Tm-T)J 

For example, the x-location of the second data point was calculated using 


x = log 


(2369 - 75 ) 


(2369 - 600 ) 


= .113 , 


where 2369° F is the melting temperature and 600° F is the x-axis value of the 2nd data 
point in Fig 4 . The y-axis values were calculated by taking the log of each yield strength 
value from Fig. 4. Referring once again to the second data point, the y-locaaon was 
calculated using y = log 156 = 2.193, where 156 ksi is the y-axis value of the 2nd data 

point in Fig. 4. 
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7.4.2 Mechanical Fatigue Data 

Figure 7 shows the effect of cycles on fatigue strength for INCONEL 718 given a 
set testing temperature. As seen from the figure, the fatigue strength of the material 
decreases as the number of cycles increases. 

Figure 8 is a plot of the data from Figure 7 that has been normalized along both 
axes. The data along the x-axis were normalized by dividing the fatigue cycle values by the 
material's fatigue cycle value at minimum fatigue strength. Thus, the x-axis values were 
calculated from 

’ Fatigue Cycles 

.Fatigue Cycles at Minimum Strength. 


For example, the x-location of the third data point was calculated using 


x 


10000 X 10 * 
100000 x 10 3 


0.1 , 


where (10000 x 10 3 ) and (100000 x 10 3 ) cycles are the x-axis values of the 3rd and 4th 
data points, respectively, in Fig. 7. Note that the x-axis value of the 4th data point in 
Fig. 7 is the above-mentioned material’s fatigue cycle value at minimum fatigue strength. 
The data along the y-axis were normalized by dividing the material’s fatigue strength values 
by its fatigue strength value at 0.1% of the fatigue cycle value at minimum fatigue strength. 
For example, the y-location of the third data point was calculated using 


where 1 1 1 and 95 ksi are the y-axis values of the 1st and 3rd data points, respectively, in 
Fig. 7 . The 1 1 1 ksi value (y-axis value of the 1st data point in Fig. 7) is the material's 
fatigue strength at 0.1% of the fatigue cycles at minimum strength, or in other words, 1 1 1 

ksi is the material's fatigue strength at 0.1% x (100000 x 10 3 ) or (100 x 10 3 ) cycles. 


Figure 9 is a plot of the log of the same data (Fig. 7) with only the x-axis 
normalized. The data along the x-axis were normalized with respect to FC m ax and FCmin, 
where FC m ax is the material's fatigue cycle value at minimum fatigue strength and FCmin 
is 0.1% of the fatigue cycle value at minimum fatigue strength. In order to calculate the 
normalized values, the difference between FC max and FCmin was computed. This value 
was then divided by the difference between FCmax and each fatigue cycle value, FC, from 
Fig. 7. Finally, the log of each normalized value was plotted along the x-axis. Thus, the 
x-axis values were calculated from 


x = log 


(FCmix ~ FCmin) 

. (FCbom - FC) . • 


For example, the x-location of the third data point was calculated using 


x = log 


(lOOOOOxlO 3 - lOOxlO 3 ) 


aoooooxio 3 - iooooxio 3 ) 


= .0453 


where (100 x 10 3 ), (10000 x 10 3 ) and (100000 x 10 3 ) are the x-axis values of the 1st, 3rd 
and 4th data points, respectively, in Fig. 7 . The y-axis values were calculated by taking 
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the log of each fatigue strength value. Referring once again to the third data point, the y- 
location was calculated using y = log 95 = 1.978 , where 95 ksi is the y-axis value of the 
3rd data point in Fig. 7.2.1. 
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Figure 8 Effect of Cycles on Fatigue Strength 
for INCONEL 718, Normalized Data 
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Figure 9 Effect of Cycles on Fatigue Strength 
for INCONEL 718, Log of Renormalized Data 



7.4.3 Creep Data 

Figure 10 shows the effect of time on rupture strength for INCONEL 718 given a 
set testing temperature. As seen from the figure, the rupture strength of the material 
decreases as the rupture life increases. 

Figure 1 1 is a plot of the data from Figure 10 that has been normalized along both 
axes. The data along the x-axis were normalized by dividing the rupture life values by the 
material's rupture life value at minimum rupture strength. Thus, the x-axis values were 
calculated from 

x Rupture Life 1 

[Rupture Life at Minimum Strength] 

For example, the x-locadon of the sixth data point was calculated using 

x* [-24210.1 = 0394 . 

L21523.6J 

where 8473.0 and 21523.6 hours are the x-axis values of the 6th and 7th data points, 
respectively, in Fig. 10 . Note that the x-axis value of the 7th data point in Fig. 7 is the 
above-mentioned material’s rupture life value at minimum rupture strength. The data along 
the y-axis were normalized by dividing the material's rupture strength values by its rupture 
strength value at 0.13% of the rupture life at minimum rupture strength. For example, the 
y-location of the sixth data point was calculated using 

y = [l2il = ,7 85 , 

L 158 J 

where 158 and 124 ksi are the y-axis values of the 1st and 6th data points, respectively, in 
Fig. 10 . The 158 ksi value (y-axis value of the 1st data point in Fig. 10) is the material's 
rupture strength at 0.13% of the rupture life at minimum strength, or in other words, 158 

ksi is the material's rupture strength at (0.13% x 21523.6) or 27.8 hours. 

Figure 12 is a plot of the log of the same data (Fig. 10) with only the x-axis 
normalized. The data along the x-axis were normalized with respect to RL ma x and RLmin, 
where RLmax is the material's rupture life value at minimum rupture strength and RL m j n is 
0.13% of the rupture life value at minimum rupture strength. In order to calculate the 
normalized values, the difference between RLmax and RLmin was computed. This value 
was then divided by the difference between RL max and each rupture life value, RL, from 
Fig. 10. Finally, the log of each normalized value was plotted along the x-axis. Thus, the 
x-axis values were calculated from 


x = log 


(RLmax * RLmin) 

- (RLmax -RL) - . 


For example, the x-location of the sixth data point was calculated using 


x = log 


(21523.6 - 27.8) 


(21523.6 - 8473.0 


= .217 , 


where 27.8, 8473.0 and 21523.6 are the x-axis values of the 1st, 6th and 7th data points, 
respectively, in Fig. 10 . The y-axis values were calculated by taking the log of each 
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rupture strength value. Referring once again to the sixth data point, the y-location was 
calculated using y = log 124 = 2.093, where 124 ksi is the y-axis value of the 6th data 
point in Fig. 10 . 
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Figure 1 0 Effect of Time on Rupture Strength 
for INCONEL 718, Raw Data 
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Figure 1 1 Effect of Time on Rupture Strength 
for INCONEL 718, Normalized Data 
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8.0 USING BOTH PROMISS AND PROMISC CODES 


8.1 Introduction 

The PROMISC code can easily be envisioned as a preprocessor to the PROMISS 
code since it provides the material constants, ai, either from expert opinion or from 
experimental data. These constants would then be input to the PROMISS code and used to 
predict strength. A number of precautions should be exercised by investigators using these 
codes in this manner. These precautions are explained in die following section. 

8.2 PROMISC - PROMISS Interface 

Before using PROMISC-gcnerated material constants as input to PROMISS, the 
diagnostics provided in the PROMISC output must be checked for violation of 
assumptions. For example, the user needs to verify that die R-squared value and the p- 
value for regression in the ANOVA Table are adequate to justify the validity of the 
regression model. If the value of R-squared is too small then the multiple linear regression 
equation used in PROMISC has little or no predictive capability. Then, the user needs to 
glean insignificant effects or primitive variables from the multifactor interaction equation. 
The significance of an effect can be determined from the column of p-values in the 
Sequential Statistics Table. If the user does not remove insignificant variables, then the 
multifactor interaction equation, used in PROMISS, would have an extra product term due 
to an extraneous regressor variable. Also, the user needs to scrutinize the Case Analysis 
Table for the presence of outliers by examining the column which contains Cook's 
Distance. The presence of an outlier can inordinately influence the regression coefficients 
which are the exponents in the multifactor interaction equation (i.e., material constants, aj). 
If these regression coefficients change markedly, then so will the PROMISS-predicted 
strength. Finally, the user should test the residuals from the Case Analysis Table to see if 
they are normally distributed. This analysis can be completed through the Shapiro-Wilk 
test If the residuals are not normal, then a transformation of the regression variables 
(effects or primitive variables) may be indicated. Once the user has positively answered the 
previous four caveats, then predicting strength using PROMISS from the multiple linear 
regression equation in PROMISC is valid. Further details are available in Section 1 1 .0, 
APPENDIX 2 (Sections 3.1.1 and 3.2.1, Discussion of Results). 
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1.0 INTRODUCTION 


This User Manual documents the FORTRAN program PROMISS (Vers. 2.0). The 
program determines the random strength of an aerospace propulsion component, due to a 
number of diverse random effects (see Section 2.0, Theoretical Background). 

Included in this User Manual are details regarding the theoretical background of 
PROMISS, input data instructions and sample problems illustrating the use of both 
PROMISS and PROMISC. APPENDIX A gives information on the primitive variables, 
their symbols, FORTRAN names and both SI and U.S. Customary units. APPENDIX B 
includes a disk containing the actual input and output files corresponding to the sample 
problems. The source code is available from the first author at the address given on the 
cover page of this report APPENDIX C details the IMSL, Version 10 [1], subroutines 
and functions called by PROMISS. APPENDIX D illustrates SAS/GRAPH [2] programs 
that can be used to plot both the probability density functions (p.d.f.) and the cumulative 
distribution functions (c.d.f.). 
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2.0 THEORETICAL BACKGROUND 


Recently, a general phenomenological constitutive relationship, for composite 
materials subjected to a number of diverse effects or primitive variables, has been 
postulated to predict mechanical and thermal material properties [3, 4, 5, 6]. The resulting 
multifactor interaction constitutive equations summarize composite micromechanics theory 
and have been used to predict material properties for a unidirectional fiber-reinforced 
lamina, based on the corresponding properties of the constituent materials. 

These equations have been modified to predict the mechanical property of strength 
for one constituent material due to "n” diverse effects or primitive variables. These effects 
could include both time-independent and time-integrated primitive variables, such as 
mechanical stresses subjected to both static and impact loads, thermal stresses due to 
temperature variations and thermal shock, and other effects such as chemical reaction or 
radiation attack. They might also include other tune-dependent primitive variables such as 
creep, mechanical fatigue, thermal aging, thermal fatigue, or even effects such as seasonal 
attack (see APPENDIX A, Primitive Variables, Symbols, and Units). For most of these 
primitive variables, strength has been observed to decrease with an increase in the variable. 

The postulated constitutive equation accounts for the degradation of strength due to 
these primitive variables. The general form of the equation is 




A iF' A » 


(i) 


where Aj, Atf? and Ajo are the current, ultimate and reference values of a particular effect, 
aj is the value of an empirical constant for the i* effect or primitive variable, n is the 
number of product terms of primitive variables in the model, and S and Sq are the current 
and reference values of material strength. Each term has the property that if the current 
value equals the ultimate value, the current strength will be zero. Also, if the current value 
equals the reference value, the term equals one and strength is not affected by that variable. 


This deterministic constitutive model may be calibrated by an appropriately curve-fitted 
least squares multiple linear regression of experimental data [7], perhaps supplemented by 
expert opinion. Ideally, experimental data giving the relationship between effects and 
strength is obtained. For example, data for just one effect could be plotted on log-log 
paper. A good fit for the data is then obtained by a linear regression analysis. This is 
illustrated schematically in Figure 1. The postulated constitutive equation, for a single 


effect, is then obtained by noting the linear relation between log S and log 


Af- Aq 
La f -a. 
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Fig. 1 Schematic of experimental data illustrating 
the effect of one primitive variable on strength. 



as follows: 



Note that the above equation (2) is for a primitive variable that lowers strength. If a 
variable raises strength, the exponent is negative. 

This general constitutive model may be used to estimate the strength of an aerospace 
propulsion system component under the influence of a number of diverse effects or 
primitive variables. The probabilistic treatment of this equation includes randomizing the 
deterministic multifactor interaction constitutive equation, performing probabilistic analysis 
by simulation and generating probability density function (p.d.f.) estimates for strength 
using a non-parametric method, maximum penalized likelihood [8,9]. Integration yields 
the cumulative distribution function (c.d.f.) from which probability statements regarding 
strength may be made. This probabilistic constitutive model predicts the random strength 
of an aerospace propulsion component due to a number of diverse random effects. 

This probabilistic constitutive model is embodied in two FORTRAN programs, 
PROMISS (Prob abilistic Material Strength Simulator) and PROMISC (Pro babilistic 
Material Strength Calibrator); see Final Technical Report, APPENDIX 2. PROMISS 
calculates the random strength of an aerospace propulsion component due to as many as 
eighteen diverse random effects. Results are presented in the form of probability density 
functions and cumulative distribution functions of normalized strength, S/So- PROMISC 
calculates the values of the empirical material constants, aj. 

PROMISS includes a relatively simple "fixed" model as well as a "flexible" model. 
The fixed model postulates a probabilistic constitutive equation that considers the primitive 
variables given in Table 1. liie general form of this constitutive equation is given in 
equation (1), wherein there are now n = 7 product terms, one for each effect or primitive 
variable listed above. Note that since this model has seven primitive variables, each 
containing four values of the variable, it has a total of twenty-eight variables. The flexible 
model postulates a probabilistic constitutive equation that considers up to as many as n = 18 
product terms for primitive variables. These variables may be selected to utilize the theory 
and experimental data currently available for the specific strength degradation mechanisms 
of interest The specific effects included in the flexible model are listed in Table 2. Note 
that in order to provide for future expansion and customization of the flexible model, six 
"other" effects have been provided. 
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Table 1 Primitive variables available in the fixed model 


itit Primitive 
Variable 

Primitive 
Variable Type 

1 

Stress due to static load 

2 

Temperature 

3 

Chemical reaction 

4 

Stress due to impact 

5 

Mechanical fatigue 

6 

Thermal fatigue 

7 

Creep 

Table 2 Primitive variables available in the flexible model 


A. Environmental Effects 

1. Mechanical 

a. Stress 

b. Impact 

c. Other Mechanical Effect 

2. Thermal 

a. Temperature Variation 

b. Thermal Shock 

c. Other Thermal Effect 

3. Other Environmental Effects 

a. Chemical Reaction 

b. Radiation Attack 

c. Other Environmental Effect 

B. Time-Dependent Effects 

1. Mechanical 

a. Creep 

b. Mechanical Fatigue 

c. Other Mech. Time-Dep. Effect 

2. Thermal 

a. Thermal Aging 

b. Thermal Fatigue 

c. Other Thermal Time-Dep. Effect 

3. Other Time-Dependent Effects 

a. Corrosion 

b. Seasonal Attack 

c. Other Time-Dep. Effect 
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The considerable scatter of experimental data and the lack of an exact description of 
the underlying physical processes for the combined mechanisms of fatigue, creep, 
temperature variations, and so on, make it natural, if not necessary to consider probabilistic 
models for a strength degradation model. Therefore, the fixed and flexible models 
corresponding to equation (1) are ’’randomized", and yield the "random normalized material 
strength due to a number of diverse random effects or primitive variables." Note that for 
the fixed model, equation (1) has the following form: 


S/So = f(AiF, Ai, Aio, A2F, A2, A20.-~.A7F, A7, A70) (3) 


where A,, Ajf and Ajo are the ultimate, current and reference values of the i* of seven 
effects or primitive variables as given in Table 1. In general, this expression can be written 
as. 


S/So = f(Xi), i = 1 28, (4) 


where the Xj are the twenty eight independent variables in equation (3). Thus, the fixed 

model is "randomized" by assuming all the independent variables, Xi, i = 1 28, to be 

random and stochastically independent. For die flexible model, equation (1) has a form 
analogous to equations (3) and (4), except that there arc as many as seventy-two 
independent variables. Applying probabilistic analysis to either of these randomized 
equations yields the distribution of the dependent random variable, normalized material 
strength, S/So* 

Although a number of methods of probabilistic analysis are available, [8] simulation 
was chosen for PROMISS. Simulation utilizes a theoretical sample generated by numerical 
techniques for each of the independent random variables. One value from each sample is 
substituted into the functional relationship, equation (3), and one realization of normalized 
strength, S/So, is calculated. This calculation is repeated for each value in the set of 
samples, yielding a distribution of different values for normalized strength. 

A probability distribution function is generated from these different values of 
normalized strength, using a non-parametric method, maximum penalized likelihood. 
Maximum penalized likelihood generates the p.d.f. estimate using the method of maximum 
likelihood together with a penalty function to smooth it [9]. Finally, integration of the 
generated p.d.f. results in the cumulative distribution function, from which probabilities of 
normalized strength can be directly observed. 

PROMISS includes computational algorithms for both the fixed and the flexible 
probabilistic constitutive models. As described above, PROMISS randomizes the 
following equation: 


s n 

f- IT 

o i=l 


j Sf‘ A i0 


(5) 
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is the i* effect, Ai, Ajf and Aio are random variables, aj is the i* empirical material 
constant and S/So is normalized strength. There arc a maximum of eighteen possible 
effects or primitive variables that may be included in the model. For the flexible model 
option, they may be chosen by the user from those in Table 2. For the fixed model option, 
the primitive variables of Table 1 are chosen. Within each primitive variable term the 
current, ultimate and reference values and the empirical material constant may be modeled 
as either deterministic (empirical, calculated by PROMISC), normal, lognormal, or Wiebull 
random variables. Simulation is used to generate a set of realizations for normalized 
random strength, S/So, from a set of realizations for primitive variables and empirical 
material constants. Maximum penalized likelihood is used to generate an estimate for the 
p.d.f. of normalized strength, from a set of realizations of normalized strength. Integration 
of the p.clf. yields the c.d.f. Plot files are produced to plot both the p.d.f. and the c.di. 
PROMISS also provides information on S/So statistics (mean, variance, standard deviation 
and coefficient of variation). A resident database, for database rather than user input of 
empirical material constants, is also provided. 

PROMISC (see Final Technical Report, APPENDIX 2) performs a multiple linear 
regression on actual experimental or simulated experimental data for as many as eighteen 
effects or primitive variables, yielding regression coefficients that are the empirical material 
constants, a;, required by PROMISS. It produces the multiple linear regression of the log 
transformation of equation (3), the PROMISS equation. When transformed it becomes 



or 

log S = log So ♦ T a log [ A f 1 , (7) 

m LAiF - AioJ 


where 


i V‘ A i 


Atf-Agj 


is the i^ effect, Ai, AjF and Ajo are primitive variable data and aj is the i* empirical 
material constant, or the i 1 * 1 regression coefficient to be predicted by PROMISC. Also, 
log So is the log transformed reference value of strength, or the intercept regression 
coefficient to be predicted by PROMISC, and log S is the log transformed strength. 



Experimental data for up to eighteen possible effects, as given in Table 2, may be included. 
The primitive variable data may be either actual experimental data or expert opinion, directly 
read from input, or simulated data where expert opinion is specified as the mean and 
standard deviation of a normal or lognormal distribution. The simulated data option for 
input data was used in the early stages of code development to verify correct performance. 
The input data, whether actual or simulated, is read in and assembled into a data matrix. 
From this data matrix, a corrected sums of squares and crossproducts matrix is computed. 
From this sums of squares and crossproducts matrix, and a least squares methodology, a 
multiple linear regression is performed to calculate estimates for the empirical material 
constant, a;, and die reference strength. So- These are the regression coefficients. 

PROMTS C includes enhancements of die multiple linear regression analysis to 
screen data from "outliers" and collinearities, determine "how well" the data fit the 
regression, quantify the importance and relative importance of each factor in the postulated 
constitutive equation, eq. (1), as well as check assumptions inherent in the use of multiple 
linear regression. Further details are provided in the Final Technical Report, Section 6.0, 
NASA Grant No. NAG 3-867, Supp. 2, "Probabilistic Lifetime Strength of Aerospace 
Materials via Computational Simulation." 
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3.0 SAMPLE PROBLEMS, DATA INPUT, AND DISCUSSION OF RESULTS 


Data input for PROMISS is user friendly and easy to enter. Data can be directly 
input by the user or the user can select data to be input by the program from its own 
resident database. Several examples follow (see also Section 6.0, APPENDIX B). 


3. 1 PROMISS Input for Flexible Model with No Statistical Distribution (Deterministic) 
and No use of the Resident Database 

The 1st line of input (format 2E12.4, see item 1, below) determines the random 
number generator seed and the data sample size. The 2nd line (format D, see item 2, 
below) determines either a fixed or a flexible model. The 3rd through 8th lines of input 
(format I3,2X,I3,2XJ3, see item 3, below) choose the 18 effects for the flexible model. 
On the 9th line of input (format 13, see item 4, below), the user chooses if the data is to be 
read from the user input or the database. The 10th line (format I3 f 2X,I3,2X,I3, 2XJ3, 
see item 4b below) determines the nature of the effect or primitive variable (deterministic or 
random). The 1 1th through 14th lines (format 10X.D12.4, see items 4b to 22, below) 
specify values for the effect or primitive variable. This sequence of five input lines (i.e., 9 
to 14) repeats for the other effects selected (see items 5 to 21). A table listing the primitive 
variables, their units and symbols is given in Section 5.0, APPENDIX A. Finally, IMSL 
subroutine parameters are entered as indicated in items 23 and 24. 


1. Line 1 selects the Random Number Generator Seed (ISEED) and Sample Size (NTOT). 
EXAMPLE: 

12345678901234567890 * 

1 40 


2. Line 2 selects either Fixed or Flexible Model ( MODEL — 0 is flexible, MODEL — 1 is 
fixed). 

EXAMPLE: 

12345678901234567890 

0 


3. Lines 3 to 8 select the 18 effects for the Flexible Model where: 

EFFMS = 1 is Quasi-static Stress Effect 
EFFMS = 0 is No Quasi-static Stress Effect 
EFFMI = 1 is Impact Effect 
EFFMI = 0 is No Impact Effect 


* NOTE: the ruler is to aid the user in formatting and is not a part of the input. 
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EFFMO = 1 is Other Mechanical Effect 
EFFMO = 0 is No Other Mechanical Effect 

EFFTT = 1 is Temperature Effect 
EFFTT = 0 is No Temperature Effect 
EFFTS — 1 is Thermal Shock Effect 
EFFTS = 0 is No Thermal Shock Effect 
EFFTO » 1 is Other Thermal Effect 
EFFTO = 0 is No Other Thermal Effect 

EFFOC = 1 is Chemical Effect 
EFFOC = 0 is No Chemical Effect 
EFFOR = 1 is Radiation Effect 
EFF OR = 0 is No Radiation Effect 
EFFOO = 1 is Other Effect 
EFFOO = 0 is No Other Effect 

TEFFMC = 1 is Creep Effect 

TEFFMC = 0 is No Creep Effect 

TEFFMF = 1 is Mechanical Fatigue Effect 

TEFFMF = 0 is No Mechanical Fatigue Effect 

TEFFMO = 1 is Other Time-Dependent Mechanical Effect 

TEFFMO = 0 is No Other Time-Dependent Mechanical Effect 

TEFFTA = 1 is Thermal Aging Effect 

TE FFTA = 0 is No Thermal Aging Effect 

T EFFTF = 1 is Thermal Fatigue Effect 

TE FFT F = 0 is No Thermal Fatigue Effect 

TE FFT O = 1 is Other Time-Dependent Thermal Effect 

TEFFTO = 0 is No Other Time-Dependent Thermal Effect 

TEFFOC = 1 is Corrosion Effect 
TEFFOC = 0 is No Corrosion Effect 
TEFFOS = 1 is Seasonal Attack Effect 
TEFFOS = 0 is No Seasonal Attack Effect 
TEFFOO = 1 is Other Time-Dependent Effect 
TEFFOO = 0 is No Other Time-Dependent Effect 

The effects arc read from the input file as follows: 

EFFMS, EFFMI, EFFMO 
EFFTT, EFFTS, EFFTO 
EFFOC, EFFOR, EFFOO 
TEFFMC, T EFFMF , TEFFMO 
TEFFTA, TEFFTF, TEFFTO 
TEFFOC, TEFFOS, TEFFOO. 
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EXAMPLE: 


12345678901234567890 

111 

111 

111 

111 

111 

111 


4a. Line 9 selects either User Input or Resident Database for the Exponent for Quasi-static 
Stress Effect (DATA = 0 is user input , DATA = 1 is database). 

EXAMPLE: 

12345678901234567890 

0 


4b. Line 10 specifies the nature of the four variables within the effect by using flags. The 
flag names are AFDS, ADS, AODS, and SADS. Setting a flag to 0 indicates a 
deterministic or nominal value will be input for the variable. Setting a flag to 1 indicates a 
normal distribution. Setting a flag to 2 indicates a lognormal distribution. Setting a flag to 
3 indicates a special distribution not yet developed. Setting a flag to 4 indicates a Weibull 
distribution. Lines 1 1 to 14 give the values of the four variables. The Quasi-static Stress 
Effect variable names have FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

1234567890123456789 01234567890 
0 0 0 0 

5.0D+00 
3 . 855D+00 
4 . 0D+00 
1.0D+00 


5. Lines 15 to 20 input data for the Impact Effect in the same manner as was described for 
the Quasi-static Stress Effect (see items 4a and 4b, above). The Impact Effect variable 
names have FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

1 23456789012345678901234567890 
0 

0 0 0 0 

5.0D+00 
3 . 855D+00 
4 . 0D+00 
1 .0D+00 
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6. Lines 21 to 26 input data for the Other Mechanical Effect in the same manner as was 
described for the Quasi-static Stress Effect The Other Mechanical Effect variable names 
have FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

1234567890123456789 01234567890 

0 

0 0 0 0 

5.0D+00 
3 . 855D+00 
4.0D+00 
1 .0D+00 


7. Lines 27 to 32 input data for the Temperature Effect in the same manner as was 
described for the Quasi-static Stress Effect. The Temperature Effect variable names have 
FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

123456789012345678901234567890 

0 

0 0 0 0 

5 . 0D+00 
3 . 855D+00 
4 . 0D+00 
1 .0D+00 


8. Lines 33 to 38 input data for the Thermal Shock Effect in the same manner as was 
described for the Quasi-static Stress Effect. The Thermal Shock Effect variable names have 
FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

123456789012345678901234567890 

0 

0 0 0 0 

5.0D+00 
3 . 855D+00 
4 . 0D+00 
1 . 0D+00 
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9. Lines 39 to 44 input data for the Other Thermal Effect in the same manner as was 
described for the Quasi-static Stress Effect The Other Thermal Effect variable names have 
FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

1 7 34567890 1 2345678901 234 5678 9.Q 
0 

0 0 0 0 

5.0D+00 
3 . 855D+00 
4 . 0D+00 
1 .0D+00 


10. Lines 45 to 50 input data for the Chemical Effect in the same manner as was described 
for the Quasi-static Stress Effect. The Chemical Effect variable names have FORTRAN 
names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

12345678901234567890 1234567890 

0 

0 0 0 0 

5 . 0D+00 
3 . 855D+00 
4 . 0D+00 
1 .0D+00 


1 1. Lines 51 to 56 input data for the Radiation Effect in the same manner as was described 
for the Quasi-static Stress Effect. The Radiation Effect variable names have FORTRAN 
names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

123456789012345678 901234567890 

0 

0 0 0 0 

5.0D+00 
3 . 855D+00 
4 . 0D+00 
1.0D+00 
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12. Lines 57 to 62 input data for the Other Effect in the same manner as was described for 
the Quasi-static Stress Effect The Other Effect variable names have FORTRAN names 
given in Section 5.0, APPENDIX A. 

EXAMPLE: 

12345 67 89012345 67 8 90 12345 67 8 9Q 
0 

0 0 0 0 

5.0D+00 
3 . 855D+00 
4 . 0D+00 
1 .0D+00 


13. Lines 63 to 68 input data for the Creep Effect in the same manner as was described for 
the Quasi-static Stress Effect. The Creep Effect variable names have FORTRAN names 
given in Section 5.0, APPENDIX A. 

EXAMPLE: 

123456789012345678901234567890 

0 

0 0 0 0 

5.0D+00 
3 . 855D+00 
4 . 0D+00 
1.0D+00 


14. Lines 69 to 74 input data for the Mechanical Fatigue Effect in the same manner as was 
described for the Quasi-static Stress Effect The Mechanical Fatigue Effect variable names 
have FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

123456789012345678901234567890 

0 

0 0 0 0 

5.0D+00 
3 . 855D+00 
4 . 0D+00 
1.0D+00 
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15. Lines 75 to 80 input data for the Time-dependent Mechanical Effect in the same 
manner as was described for the Quasi-static Stress Effect. The Tune-dependent 
Mechanical Effect variable names have FORTRAN names given in Section 5.0, 
APPENDIX A. 

EXAMPLE: 

1 734.S6789Q1 7345 678 90 1 934567890 
0 

0 0 0 0 

5.0D+00 
3 . 855D+00 
4 . 0D+00 
1.0D+00 


16. Lines 81 to 86 input data for the Thermal Aging Effect in the same manner as was 
described for the Quasi-static Stress Effect The Thermal Aging Effect variable names have 
FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

1 ?34 5 6789012345 67 8901234567890 
0 

0 0 0 0 

5 . 0D+00 
3 . 855D+00 
4 . 0D+00 
1.0D+00 


17. Lines 87 to 92 input data for the Thermal Fatigue Effect in the same manner as was 
described for the Quasi-static Stress Effect The Thermal Fatigue Effect variable names 
have FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

1 234 5 678901234 5 67 89012345678-SO 
0 

0 0 0 0 

5 . 0D+00 
3 . 855D+00 
4 . 0D+00 
1 .0D+00 
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18. Lines 93 to 98 input data for the Other Time-dependent Thermal Effect in the same 
manner as was described for the Quasi-static Stress Effect The Other Time-dependent 
Thermal Effect variable names have FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

123456789012345678901234567890 

0 

0 0 0 0 

5.0D+00 
3 . 855D+00 
4 . 0D+00 
1 .0D+00 


19. Lines 99 to 104 input data for the Corrosion Effect in the same manner as was 
described for the Quasi-static Stress Effect The Corrosion Effect variable names have 
FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

123456789012345678901234567890 

0 

0 0 0 0 

5.0D+00 
3 . 855D+00 
4 . 0D+00 
1 .0D+00 


20. Lines 105 to 1 10 input data for the Seasonal Attack Effect in the same manner as was 
described for the Quasi-static Stress Effect The Seasonal Attack Effect variable names 
have FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

123456789012345678901234567890 

0 

0 0 0 0 

5.0D+00 
3 . 855D+00 
4 . 0D+00 
1 .0D+00 
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2 1 Lines 1 1 1 to 1 1 6 input data for the Other Time-dependent Effect in the same manner as 
was described for the Quasi-static Stress Effect. The Other Time-dependent Effect variable 
names have FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

1? 3 4567890 1734 5 678 90 12345678 9Q 
0 

0 0 0 0 

5.0D+00 
3 . 855D+00 
4 . 0D+00 
1.0D+00 


22. The DESPL [1] parameters are NODE, INTT, ALPHA, EPS, and MAXIT and are 
entered in that order as follows (line 1 17): 

EXAMPLE: 

1234567890123456789012345678901234567890 
21 0 1.0E+01 1.0E-05 30 

23. The DESPL[ 1] parameter, IOPT, is entered as follows (line 118): 

EXAMPLE: 

1234567890 

2 


3.1.1 Discussion of Results 

Execution of PROMISS (source code entitled PROMISS89.FOR ) produces an 
output file that gives numerical results (see Section 6.0, APPENDIX B). Execution also 
produces plotfiles (see Section 6.0, APPENDIX B). These files are used to plot the X and 
Y axes of the probability density function (p.d.f.) and the cumulative distribution function 
(c.d.f.) generated by PROMISS. The plots are drawn from the plotfiles by the 
SAS/GRAPH graphing program (see Section 8.0, APPENDIX D). 
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3.2 PROMISS Input for Fixed Model with Various Statistical Distributions and use of 
the Resident Database 

The 1st line of input (format 2E12.4, see item 1, below) determines the random 
number generator seed and the data sample size. The 2nd line (format 13, see item 2, 
below) determines either a fixed or a flexible model On the 3rd line of input the user 
chooses if the data is to be read from the user input or die resident database (format 13, see 
item 3, below). The 4th line (format 13,2X43,2X43, 2X43, see item 4, below) 
determines the nature of the effect or primitive variable (deterministic or random). The 5th 
through 8th lines specify values for the effect or primitive variable. This sequence of five 
input lines repeats for die other effects (see items 5 through 10, below). If the word 
DATABASE is entered in the place of the fourth integer (format I3,2X43,2X43,2X,A8), 
the program will read a deterministic value for the exponent from its own Resident 
Database. A table listing the primitive variables, their units and symbols is given in 
APPENDIX A. Finally, IMSL subroutine parameters are entered as indicated in items 12 
and 13. 

1. Line 1 selects Random Number Generator Seed (ISEED) and Sample Size (NTOT). 
EXAMPLE: 

12345678901234567890 * 

1 40 


2. Line 2 determines fixed or flexible model (MODEL = 0 is flexible, MODEL = 1 is 
fixed). 

EXAMPLE: 

12345678901234567890 

1 


3. Line 3 determines User Input or Database for the Exponent for Quasi-Static Stress 
Effect (DATA = 0 is user input, DATA = 1 is database). 

EXAMPLE: 

12345678901234567890 

1 


* NOTE: the ruler is to aid the user in formatting and is not a part of the input. 
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4. Line 4 specifies the nature of the four variables within the effect by using flags. The flag 
names are AFDS, ADS, AODS, and SADS. Setting a flag to 0 indicates a deterministic or 
nominal value will be input for the variable. Setting a flag to 1 indicates a normal 
distribution. Setting a flag to 2 indicates a lognormal distribution. Setting a flag to 3 
indicates a special distribution not yet developed. Setting a flag to 4 indicates a Weibull 
distribution. Lines 5 to 7 give the values of die four variables. The Quasi-static Stress 
Effect variable names have FORTRAN names given in Section 5.0, APPEND DC A. 

EXAMPLE: 

1234567890123456789012345678901234567890 
4 4 2 DATABASE “ ” 

130.0D+00 6 . 500D+00 

90.0D+00 4 . 500D+00 

-2.9D+00 -0 . 145D+00 


5. Line 8 determines user input or database for the exponent for Impact Effect Lines 9 to 
13 input data for the Impact effect in the same manner as was described for the Quasi-static 
Stress Effect (see items 3 and 4, above). The Impact Effect variable names have 
FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

1234567890123456789012345678 901234567890 

0 

4 12 4 

1.0D+00 0 . 003D+00 

0.10D+00 0 . 003D+00 

0 . 001D+00 0 . 00003D+00 
-0.5D+00 -0 . 015D+00 


6. Line 14 determines user input or database for the exponent for Temperature Effect. 
Lines 15 to 19 input data for the Temperature Effect in the same manner as was described 
for the Quasi-static Stress Effect. The Temperature Effect variable names have FORTRAN 
names given in Section 5.0, APPENDIX A. 

EXAMPLE: 


12445678901234567890123456789 01234567890 


0 

2 


1 2 
2732 . 0D+00 
1562 . 0D+00 
68.0D+00 
0.50D+00 


82 . 000D+00 
46.70D+00 
2 . 040D+00 
0 . 015D+00 
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7. Line 20 determines user input or database for die exponent for Chemical Reaction 
Effect Lines 2 1 to 25 input data for the Chemical Reaction Effect in the same manner as 
was described for the Quasi-static Stress Effect The Chemical Reaction Effect variable 
names have FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 


1234567890123456789 012345678901234567890 

0 

1111 

1.0D+00 0 . 003D+00 

0.02D+00 0 . 0006D+00 

0 . 001D+00 0 . 00003D+00 
0.50D+00 0 . 015D+00 


8. Line 26 determines user input or database for die exponent for Creep Effect Lines 27 
to 30 input data for the Creep Effect in the same manner as was described for the Quasi- 
static Stress Effect The Creep Effect variable names have FORTRAN names given in 
Section 5.0, APPENDIX A. Note that when the DATABASE option is used, the input is 
only 5 lines long per effect rather than 6 lines (compare to item 7, above). 


EXAMPLE: 

1234567890123456789012345678901234567890 

1 

242 DATABASE 

10000. 0D+00 500 . 00D+00 

105.0D+00 3.15D+00 

0 . 083D+00 0 . 0025D+00 


9. Line 3 1 determines user input or database for the exponent for Mechanical Fatigue 
Effect Lines 32 to 36 input data for the Mechanical Fatigue Effect in the same manner as 
was described for the Quasi-static Stress Effect The Mechanical Fatigue Effect variable 
names have FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 


1234567890123456789012345678901234567890 


0 

2 


2 1 

7 . 0D+00 
3.5D+00 
1 .0D+00 
0.50D+00 


0.700D+00 
0 .350D+00 
0 . 100D+00 
0 .015D+00 
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10. Line 37 determines user input or database for the exponent for Thermal Fatigue Effect. 
Lines 38 to 41 input data for the Thermal Fatigue Effect in the same manner as was 
described for the Quasi-static Stress Effect The Thermal Fatigue Effect variable names 
have FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

12345 67 89012345 67 89012345678 9Q 12345 67 89Q 
1 

222 DATABASE 

3.0D+00 0 . 300D+00 

2.3D+00 0 . 230D+00 

1.0D+00 0 . 100D+00 


11. The DESPL[ 1 ] parameters are NODE, INIT, ALPHA, EPS, and MAXIT and are 
entered in that order as follows (line 42): 

EXAMPLE: 

1234567890123456789012345678901234501234567890 
21 0 1.0E+01 1.0E-05 30 


12. The DESPL[1] parameter, IOPT, is entered as follows (line 43): 

EXAMPLE: 

1234567890 

2 


DATABASE INPUT FILE** 

12345678901234567890 
0 . 5D+00 
0.5D+00 
0.5D+00 


3.2.1 Discussion of Results 

Execution of PROMISS (source code entitled PROMISS89.FOR ) produces an 
output file that gives numerical results (see Section 6.0, APPENDIX B). Execution also 
produces plotfiles (see Section 6.0, APPENDIX B). These files are used to plot the X and 
Y axes of the probability density function (p.d.f.) and the cumulative distribution function 
(c.d.f.) generated by PROMISS. The plots arc drawn from the plotfiles by the 
SAS/GRAPH graphing program (see Section 8.0, APPENDIX D). These plots for the 
sample problem are shown in Figures 2 and 3. 


** PROMISS, Vers. 2.0, currently has a value of 0.5 for all emperical material constants. This is cun-ent 
expert opinion. These values will be updated as data become available (see Section 6.0). 
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3.3 PROMISS Input for Flexible Model with Various Statistical Distributions and Use 
of the Resident Database 

The 1st line of input (format 2E12.4, see item 1, below) determines the random 
number generator seed and the data sample size. The 2nd line (format 13, see item 2, 
below) determines either a fixed or a flexible model. The 3rd through 8th lines (format 
I3,2X,I3,2X,I3, see item 3, below) choose the 18 effects for the flexible model. The 9th 
line (format 13, see item 4, below) the user chooses if the data is to be read from the user 
input or the database. The 10th line (format B,2X,I3,2X,I3, 2X,I3, see item 4, below) 
determines the nature of the effect or primitive variable. If the word DATABASE is entered 
in the place of the fourth integer (format I3,2XJ3,2X,I3,2X,A8), the program will read a 
deterministic value from its own Resident Database. The 1 1th through 14th lines (format 
10X.2D12.4, see item 4, below) specify values for the effect or primitive variable. This 
sequence of five input lines repeats for the other effects selected (see items 4 through 21, 
below). A table listing the primitive variables, their units and symbols is given in Section 
5.0, APPENDIX A. Finally, IMSL subroutine parameters are entered as indicated in 
items 23 and 24. 

1. Line 1 selects the Random Number Generator Seed (ISEED) and Sample Size (NTOT). 
EXAMPLE: 

12345678901234567890 * 

1 40 


2. Line 2 determines fixed or flexible model (MODEL = 0 is flexible, MODEL = 1 is 
fixed). 

EXAMPLE: 

12345678901234567890 

0 


3. Lines 3 to 8 select the 18 Effects for the Flexible Model where: 

EFFMS = 1 is Quasi-static Stress Effect 
EFFMS * 0 is No Quasi-static Stress Effect 
EFFMI = 1 is Impact Effect 
EFFMI = 0 is No Impact Effect 
EFFMO = 1 is Other Mechanical Effect 
EFFMO = 0 is No Other Mechanical Effect 

EFFTT = 1 is Temperature Effect 
EFFT T = 0 is No Temperature Effect 
EFFTS = 1 is Thermal Shock Effect 
EFFT S = 0 is No Thermal Shock Effect 
EFFTO = 1 is Other Thermal Effect 
EFFTO = 0 is No Other Thermal Effect 


* NOTE: the ruler is to aid the user in formatting and is not a part of the input 
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EFFOC = 1 is Chemical Effect 
EFFOC = 0 is No Chemical Effect 
EFFOR = 1 is Radiation Effect 
EFFOR = 0 is No Radiation Effect 
EFFOO = 1 is Other Effect 
EFFOO = 0 is No Other Effect 

TEFFMC = 1 is Creep Effect 
TEFFMC = 0 is No Creep Effect _ 

TEFFMF = 1 is Mechanical Fatigue Effect^ 

TEFFMF = 0 is No Mechanical Fatigue Effect 
TEFFMO = 1 is Other Tune-Dependent Mechanical Effect 
TEFFMO = 0 is No Other Time-Dependent Mechanical Effect 

TEFFTA = 1 is Thermal Aging Effect 
TEFFTA = 0 is No Thermal Aging Effect 
TEFFIF = 1 is Thermal Fatigue Effect 
TEFFTF = 0 is No Thermal Fatigue Effect _ 

TEFFTO = 1 is Other Time-Dependent Thermal Effect 
TEFFTO = 0 is No Other Time-Dependent Thermal Effect 

TEFFOC = 1 is Corrosion Effect 
TEFFOC = 0 is No Corrosion Effect 
TEFFOS = 1 is Seasonal Attack Effect 
TEFFOS = 0 is No Seasonal Attack Effect 
TEFFOO = 1 is Other Time-Dependent Effect 
TEFFOO = 0 is No Other Time-Dependent Effect. 

The effects are read from the input file as follows: 

EFFMS, EFFMI, EFFMO 
EFFTT, EFFTS, EFFTO 
EFFOC, EFFOR, EFFOO 
TEFFMC, TEFFMF, TEFFMO 
TEFFTA, TEFFTF, TEFFTO 
TEFFOC, TEFFOS, TEFFOO. 

EXAMPLE: 

12345678901234567890 

111 

111 

111 

111 

111 

111 
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4. Line 9 selects either User Input or Resident Database for the Exponent for Quasi-static 
Stress Effect (DATA = 0 is user input , DATA = 1 is database). Line 10 specifies the 
nature of the four variables within the effect by using flags. The flag names are AFDS, 
ADS, AODS, and SADS. Setting a flag to 0 indicates a deterministic or nominal value will 
be input for the variable. Setting a flag to 1 indicates a normal distribution. Setting a flag 
to 2 indicates a lognormal distribution. Setting a flag to 3 indicates a special distribution 
not yet developed. Setting a flag to 4 indicates a WeibuU distribution. Lines 1 1 to 14 give 
the values of the four variables. The Quasi-static Stress Effect variables names have 
FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 


1234567890123456789012345678901234567890 


0 

2 4 


4 1 

5.0D+00 
3 . 855D+00 
4 . 0D+00 
0.51D+00 


0 . 002D+00 
0 . 005D+00 
0 . 007D+00 
0 . 001D+00 


5. line 15 determines user input or database for the exponent for Impact Effect lines 16 
to 19 input data for the Impact effect in the same manner as was described for the Quasi- 
static Stress Effect (see items 3 and 4, above). The Impact Effect variable names have 
FORTRAN names given in Section 5.0, APPENDIX A. Note that when the DATABASE 
option is used, the input is 5 lines long per effect, rather than 6 lines (compare to item 4 
above). 

EXAMPLE: 

1234567890123456789012345678901234567890 

1 

412 DATABASE 

5.0D+00 0 . 002D+00 

3 . 855D+00 0 . 005D+00 

4.0D+00 0 . 007D+00 


6. Line 20 determines user input or database for the exponent for Mechanical Effect lines 
21 to 24 input data for the Mechanical Effect in the same manner as was described for the 
Quasi-static Stress Effect The Mechanical Effect variable names have FORTRAN names 
given in Section 5.0, APPENDIX A. 

EXAMPLE: 

1234567890123456789012345678901234567890 

1 

442 DATABASE 

5.0D+00 0 . 002D+00 

3 . 855D+00 0 . 005D+00 

4.0D+00 0 . 007D+00 
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7. Line 25 determines user input or database for the exponent for Temperature Effect. 
Lines 26 to 29 input data for the Temperature Effect in the same manner as was described 
for the Quasi-static Stress Effect. The Temperature Effect variable names have FORTRAN 
names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

1 7 34' 5678901 2345 6 7 ft 901 234 5678901234567830 
1 

114 DATABASE 

5.0D+00 0 . 002D+00 

3 . 855D+00 0 . 005D+00 

4.0D+00 0 . 007D+00 


8 Line 30 dete rmin es user input or database for the exponent for Thermal Shock Effect, 
lines 31 to 35 input data for the Thermal Shock Effect Effect in the same manner as was 
described for the Quasi-static Stress Effect. The Thermal Shock Effect variable names have 
FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 


1234 5 6789012345 67 890123456789012345 67 8 2.Q 


0 

4 2 


1 2 

5.0D+00 
3 . 855D+00 
4 . 0D+00 
0.51D+00 


0 .002D+00 
0 .005D+00 
0 . 007D+00 
0 . 001D+00 


9. Line 36 determines user input or database for the exponent for Other Thermal Effect 
Lines 37 to 40 input data for the Other Thermal Effect in the same manner as was described 
for the Quasi-static Stress Effect The Other Thermal Effect variable names have 
FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

1234567890123456789012345678 901234567890 

1 

421 DATABASE 

5.0D+00 0 . 002D+00 

3 . 855D+00 0 . 005D+00 

4.0D+00 0 . 007D+00 
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10. Line 4 1 determines user input or database for the exponent for Chemical Reacdon 
Effect Lines 42 to 45 input data for the Chemical Reaction Effect in the same manner as 
was described for the Quasi-static Stress Effect The Chemical Reaction Effect variable 
names have FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

1234567890 1234567890 12345678901234567890 

1 

414 DATABASE 

5.0D+00 0 . 002D+00 

3 . 855D+00 0 . 005D+00 

4.0D+00 0.007D+00 


11. Line 46 determines user input or database for the exponent for Radiation Effect. Lines 
47 to 50 input data for the Radiation Effect in the same manner as was described for the 
Quasi-static Stress Effect The Radiation Effect variable names have FORTRAN names 
given in Section 5.0, APPENDIX A. 

EXAMPLE: 

1234567890123456789012345678901234567890 

1 

414 DATABASE 

5.0D+00 0 . 002D+00 

3 . 855D+00 0 . 005D+00 

4.0D+00 0 . 007D+00 


12. Line 5 1 determines user input or database for the exponent for Other Effect. Lines 52 
to 55 input data for the Other Effect in the same manner as was described for the Quasi- 
static Stress Effect. The Other Effect variable names have FORTRAN names given in 
Section 5.0, APPENDIX A. 

EXAMPLE: 

1234567890123456789012345678901234567890 

1 

142 DATABASE 

5.0D+00 0 . 002D+00 

3 . 855D+00 0 . 005D+00 

4.0D+00 0 . 007D+00 
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13. Line 56 determines user input or database for the exponent for Creep Effect. Lines 57 
to 60 input data for the Creep Effect in the same manner as was described for the Quasi- 
static Stress Effect The Creep Effect variable names have FORTRAN names given in 
Section 5.0, APPENDIX A. The last line for this Effect is eliminated when the database is 

used. 

EXAMPLE: 

i 2 34 5678901 234 56 78 901 2345678 9Q1234 5678 9.Q 
1 

122 DATABASE 

5.0D+00 0 . 002D+00 

3 . 855D+00 0 . 005D+00 

4.0D+00 0 . 007D+00 


14. Line 61 deter min es user input or database for the exponent for Mechanical Fatigue 
Effect Lines 62 to 66 input data for the Mechanical Fatigue Effect in the same manner as 
was described for the Quasi-static Stress Effect The Mechanical Fatigue Effect variable 
names have FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 


1234 5 67 8901234567890123456789012345 67 89 Q 


0 

1 


12 1 

5.0D+00 
3 . 885D+00 
4 . 0D+00 
0.51D+00 


0 .002D+00 
0 . 005D+00 
0 .007D+00 
0 . 001D+00 


15. Line 67 determines user input or database for the exponent for Other Time-dependent 
Mechanical Effect. Lines 68 to 7 1 input data for the Other Time-dependent Mechanical 
Effect in the same manner as was described for the Quasi-static Stress Effect The Other 
Time-dependent Mechanical Effect variable names have FORTRAN names given in Section 
5.0, APPENDIX A. 

EXAMPLE: 

12345678901234567890123456 78901234567890 

1 

222 DATABASE 

5.0D+00 0 . 002D+00 

3 . 855D+00 0 . 005D+00 

4.0D+00 0 . 007D+00 
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16. Line 72 determines user input or database for the exponent for Thermal Aging Effect. 
Lines 73 to 76 input data for the Thermal Aging Effect in the same manner as was 
described for the Quasi-static Stress Effect. The Thermal Aging Effect variable names have 
FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

12345678901234567890 123456789 01234567890 

1 

111 DATABASE 

5.0D+00 0 . 002D+00 

3 . 855D+00 0 . 005D+00 

4.0D+00 0 . 007D+00 


17. line 77 determines user input or database for the exponent for Thermal Fatigue Effect. 
Lines 78 to 81 input data for the Thermal Fatigue Effect in die same manner as was 
described for the Quasi-static Stress Effect. The Thermal Fatigue Effect variable names 
have FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 


1234567890123456789012345678901234567890 


1 

4 4 


4 DATABASE 
5 . 0D+00 
3 . 855D+00 
4 . 0D+00 


0 .002D+00 
0 . 005D+00 
0 . 007D+00 


18. Line 82 determines user input or database for the exponent for Other Time-dependent 
Thermal Effect Lines 83 to 86 input data for the Other Time-dependent Thermal Effect in 
the same manner as was described for the Quasi-static Stress Effect The Other Time- 
dependent Thermal Effect variable names have FORTRAN names given in Section 5.0, 
APPENDIX A. The last line for this Effect is eliminated when the database is used. 

EXAMPLE: 


1234567890123456789012345678901234567890 


1 

4 2 


2 DATABASE 
5.0D+00 
3 . 855D+00 
4 . 0D+00 


0 . 002D+00 
0 . 005D+00 
0 . 007D+00 
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19. Line 87 determines user input or database for the exponent for Corrosion Effect 
Lines 88 to 9 1 input data for the Corrosion Effect in the same manner as was described for 
the Quasi-s tati c Stress Effect The Corrosion Effect variable names have FORTRAN 
names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

1234 5 678901 7 34 5 678901234567890123456789.0 
1 

114 DATABASE 

5.0D+00 0 . 002D+00 

3 . 855D+00 0 . 005D+00 

4.0D+00 0 . 007D+00 


20. Line 92 determines user input or database for the exponent for Seasonal Attack Effect 
Lines 93 to 96 input data for the Seasonal Attack Effect in the same manner as was 
described for the Ouasi-static Stress Effect The Seasonal Attack Effect variable names 
have FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

1 234 5 678 9012345 678901234 56789012345678-911 
1 

244 DATABASE 

5.0D+00 0 . 002D+00 

3 . 855D+00 0 . 005D+00 

4.0D+00 0 . 007D+00 


21. Line 97 determines user input or database for the exponent for Other Time-dependent 
Effect Lines 98 to 101 input data for the Other Time-dependent Effect in the same manner 
as was described for the Quasi-static Stress Effect The Other Time-dependent Effect 
variable names have FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

123 4 5 678 90 1234 5 67 8 9 0 1234 5 678 9012 34 5 67 8-9.Q 
1 

424 DATABASE 

5.0D+00 0 . 002D+00 

3 . 855D+00 0 . 005D+00 

4.0D+00 0 . 007D+00 
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22. The DESPL[1] parameters are NODE, IN IT, ALPHA, EPS, and MAXIT and are 
entered in that order as follows (line 102): 

EXAMPLE: 


12345 67 89012345 67 8 90 1234567B 90 12 345 67 8 90 
21 0 1.0E+01 1.00E-05 30 


23. The DESPL[1] parameter, IOPT, is entered as follows (line 103): 

EXAMPLE: 

1234567890 

2 


DATABASE INPUT FILE 

12345678901234567890 
0 . 5D+00 
0.5D+00 
0.5D+00 
0.5D+00 
0.5D+00 
0 . 5D+00 
0 . 5D+00 
0.5D+00 
0.5D+00 
0.5D+00 
0.5D+00 
0.5D+00 
0.5D+00 
0.5D+00 
0.5D+00 
0.5D+00 
0.5D+00 
0 . 5D+00 


3.3.1 Discussion of Results 

Execution of PROMISS (source code entitled PROMISS89.FOR ) produces an 
output file that gives numerical results (see Section 6.0, APPENDIX B). Execution also 
produces plotfilcs (see Section 6.0, APPENDIX B). These files are used to plot the X and 
Y axes of the probability density function (p.d.f.) and the cumulative distribution function 
(c.d.f.) generated by PROMISS. The plots are drawn from the plotfiles by the 
SAS/GRAPH graphing program (see Section 8.0, APPENDIX D). 
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3.4 General PROMISS Program Notes 


1 . Normalized strength, S/SO, is assured to be non-negative; any negative value 
calculated from input is set arbitrarily to zero. 


2 . 


Any effect , 


Af - A 


Ap - Aq. 


aj 


not included in the model is set to one (1). 


3 . The IMSL, Vers. 10 subroutine DESPL requires that IMSL subroutine D3SPL 
be appended to PROMISS for proper operation. 

4. Fatigue cycles should be input as log cycles rather than cycles. This assures fatigue 
affects strength calculations by yielding a fraction significantly below zero. If input 
for fatigue was in cycles, the value of the effect approaches one, thereby not affecting 
strength calculations. 
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5.0 APPENDIX A 

PRIMITIVE VARIABLES, SYMBOLS, AND UNITS 
Table A 1.2 Primitive variables, symbols, and units for PROMISS 


Primitive Variables 

Theory 

FORTRAN 

Units 

U.S 

(Effect) 

Symbol 

Name 

SI 

QUASI-STATIC 
STRESS EFFECT 





Ultimate value 

Ssf 

AF 

MPa 

ksi 

Current value 

os 

A 

MPa 

ksi 

Reference value 

OSO 

AO 

MPa 

ksi 

Material constant 

p 

SA 

Dimensionless 


IMPACT EFFECT 





Ultimate value 

Sdf 

AF2 

Dimensionless 


Current value 

od 

A2 

Dimensionless 


Reference value 

OD0 

A02 

Dimensionless 


Material constant 

s 

SA2 

Dimensionless 


OTHER 

MECHANICAL 

EFFECTS 





Ultimate value 

A3F 

AF3 

Dimensionless 


Current value 

a 3 

A3 

Dimensionless 


Reference value 

A 30 

A03 

Dimensionless 


Material constant 

a 3 

SA3 

Dimensionless 


TEMPERATURE 





Ultimate value 

Tf 

AF4 

°C 

op 

Current value 

T 

A4 

°C 

op 

Reference value 

To 

A04 

°C 

°F 

Material constant 

q 

SA4 

Dimensionless 


THERMAL 

SHOCK 





Ultimate value 

A 5 F 

AF5 

Dimensionless 


Current value 

A 5 

A5 

Dimensionless 


Reference value 

A 50 

A05 

Dimensionless 


Material constant 

as 

SA5 

Dimensionless 
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Primitive Variables Theory 
(Effect) Symbol 


FORTRAN 

Name 


SI 


Units 


U.S. 


OTHER 

THERMAL 

EFFECTS 

JUfimflte value 

AfiF 

A6 

AF6 

Dimensionless 

Dimensionless 


Current value 

A6 


Reference value 

AfiO 

A06 

Dimensionless 


Material constant 

36 

SA6 

Dimensionless 


CHEMICAL 

EFFECTS 

Ultimate value 

Rp 

AF7 

Dimensionless 


Current value 

R 

A7 

Dimensionless 


Reference value 

Ro 

A07 

Dimensionless 


Material constant 

r 

SA7 

Dimensionless 


radiation 

EFFECTS 
Ultimate value 

A8F 

AF8 

Dimensionless 


Current value 

A 8 

A8 

Dimensionless 


Reference value 

Ago 

A08 

Dimensionless 


Material constant 

38 

SA8 

Dimensionless 


OTHER EFFECTS 

Ultimate value 

A9F 

AF9 

Dimensionless 


Current value 

a 9 

A9 

Dimensionless 


Reference value 

A 90 

A09 

Dimensionless 


Material constant 

39 

SA9 

Dimensionless 


CREEP 

Ultimate value 

k:F 

AF10 

Hours 

Hours 

Current value 

tc 

A10 

Hours 

Hours 

Reference value 

kx> 

AOIO 

Hours 

Hours 

Material constant 

V 

SA10 

Dimensionless 


MECHANICAL 

FATIGUE 

Ultimate value 

Nmf 

AF11 

log of cycles 


Current value 

Nm 

All 

log of cycles 


Reference value 

NmO 

AOll 

log of cycles 


Material constant 

t 

SA11 

Dimensionless 
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Primitive Variables Theory 
(Effect) Symbol 


FORTRAN 

Name 


SI 


Units 


U.S. 


OTHER TIME 
DEPENDENT 
MECHANICAL 




Ultimate value 

A 12 F 

AF12 

Dimensionless 

Current value 

A 12 

A12 

Dimensionless 

Reference value 

Ai20 

A012 

Dimensionless 

Material constant 

ai2 

SA12 

Dimensionless 

THERMAL AGING 




Ultimate value 

Ai3F 

AF13 

Dimensionless 

Current value 

A 13 

A13 

Dimensionless 

Reference value 

A 130 

A013 

Dimensionless 

Material constant 

ai3 

SA13 

Dimensionless 

THERMAL 




FATIGUE 




Ultimate value 

Ntf 

AF14 

log of cycles 

Current value 

Nj 

A14 

log of cycles 

Reference value 

Nto 

A014 

log of cycles 

Material constant 

u 

SA14 

Dimensionless 

OTHER TIME- 




DEPENDENT 




THERMAL 




EFFECTS 




Ultimate value 

Al5F 

AF15 

Dimensionless 

Current value 

Al5 

A15 

Dimensionless 

Reference value 

A 150 

A015 

Dimensionless 

Material constant 

ai5 

SA15 

Dimensionless 

CORROSION 




Ultimate value 

Ai6F 

AF16 

Dimensionless 

Current value 

Ai6 

A16 

Dimensionless 

Reference value 

Ai60 

A016 

Dimensionless 

Material constant 

ai6 

SA16 

Dimensionless 
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Primitive Variables Theory 
(Effect) Symbol 


FORTRAN 

Name 


SI 


Units 


U.S. 


SEASO NA L 

ATTACK 


Ultimate value 

Al7F 

AF17 

Dimensionless 

Current value 

Al7 

A17 

Dimensionless 

Reference value 

a 170 

A017 

Dimensionless 

Material constant 

ai7 

SA17 

Dimensionless 


OTHER TIME- 

P E P ENP ENI 




Ultimate value 

Ai8F 

AF18 

Dimensionless 

Current value 

Ai8 

A18 

Dimensionless 

Reference value 

A 180 

A018 

Dimensionless 

Material constant 

ai8 

SA18 

Dimensionless 
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6.0 APPENDIX B 

PROMISS SAMPLE PROBLEM: INPUT AND OUTPUT FILES 


S am ple problems are discussed in sections 3.1, 3.2, and 3.3. One sample problem 
corresponds to each section. The input and output file names for each sample problem are 
listed below. The enclosed disk also includes the same files. 

3. 1 PROMISS INPUT FOR FLEXIBLE MODEL WITH NO STATISTICAL 
DISTRIBUTION (DETERMINISTIC) AND NO USE OF THE RESIDENT DATABASE 

Input File(s): 31PR.INP 
Output File: 31PR.OUT 

Plot Files: 31PR1.PLT 

31PR2.PLT 

3.2 PROMISS INPUT FOR FIXED MODEL WITH VARIOUS STATISTICAL 
DISTRIBUTIONS AND USE OF THE RESIDENT DATABASE 

Input File(s): 32PR1.INP 
32PR2.INP 
Output File: 32PR.OUT 

Plot Files: 32PR1.PLT 

32PR2.PLT 

3.3 PROMISS INPUT FOR FLEXIBLE MODEL WITH VARIOUS STATISTICAL 
DISTRIBUTIONS AND USE OF THE RESIDENT DATABASE 

Input File(s): 33PR1.INP 
33PR2.INP 
Output File: 33PR.OUT 

Plot Files: 33PR1.PLT 

33PR2.PLT 
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7.0 APPENDIX C 


PROMISS 

1. RNSET 

2. RNNOR 

3. RNLNL 

4. RNWEB 

5. DESPL 

6. GCDF 


IMSL SUBROUTINE CALLS FROM PROMISS 


Initializes a random seed for use in the IMSL 
random number generators. 

Generates pseudorandom numbers from a standard 
normal distribution using an inverse CDF 
method. 

Generates pseudorandom numbers from a 
lognormal distribution. 

Generates pseudorandom numbers from a 
Weibull distribution. 

Performs nonparametric probability density 
function estimation by the penalized 
likelihood method. 

Evaluates a general continuous cumulative 
distribution function given the ordinates 
of the density. 


7. Other IMSL Subroutines 


SSCAL and SADD 


8.0 APPENDIX D 

SAMPLE SAS/GRAPH PROGRAM FOR PROMISS 


data a; 

INFILE 'RAND1.CPR' FIRSTOBS=2;input x y; 

GOPTIONS DEVICE=HP7 470; 
proc gplot; 

axisl labcl=(h=l f=simplcx 'LOG OF CYCLES ) 
value=(h=l f=simplcx); 
axis2 valuc=(h=l f=simplcx) labcl=none; 

plot y*x / haxis=axis 1 vaxis=axis2; 

TITLE H=1 A =90 F=SIMPLEX 'PROBABILITY DENSITY FUNCTION'; 
symbol i=spline v=square; 
dataB; 

INFILE 'RAND2.CPR' FIRSTOBS=2;input x y; 
proc gplot; 

axisl label=(h=l f=simplcx 'LOG OF CYCLES') 
value=(h=l f=simplex); 
axis2 valuc=(h=l f=simplex) labcl=none; 

plot y*x / haxis=axis 1 vaxis=axis2; 

TITLE H=1 A =90 F=SIMPLEX 'CUMULATIVE DISTRIBUTION FUNCTION ; 
symbol i=spline v=square; 
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1.0 INTRODUCTION 


This User Manual documents the FORTRAN program PROMISC. The program 
performs a multiple linear regression on actual experimental or simulated experimental data 
for as many as eighteen effects or primitive variables, yielding regression coefficients that 
are the empirical material constants, ai, required by PROMISS (see Section 2.0, Theoretical 

Background). 


Included in this User Manual are details regarding the theoretical background of 
both PROMISS and PROMISC, input data instructions and sample problems illustrating 
the use of PROMISC. Appendix A gives information on the primitive variables, their 
symbols FORTRAN names and both SI and U.S. Customary units. Appendix B includes 
a disk containing the actual input and output files corresponding to the sample problems. 
The source code is available from the author at the address given on the cover page of tins 
report. Appendix C details the IMSL, Version 10 [1], subroutines and functions called by 

PROMISC. 
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2.0 THEORETICAL BACKGROUND 


Recently, a general phenomenological constitutive relationship for composite 
materials subjected to a number of diverse effects or primitive variables has been postulated 
to predict mechanical and thermal material properties [3, 4,5, 6]. The resulting multifactor 
interaction constitutive equations summarize composite micromechanics theory and have 
been used to predict material properties for a unidirectional fiber-reinforced lamina, based 
on the corresponding properties of the constituent materials. 

These equations have been modified to predict the mechanical property of strength 
for one constituent material due to "n" diverse effects or primitive variables. These effects 
could include both time-independent and time-integrated primitive variables, such as 
mechanical stresses subjected to both static and impact loads, thermal stresses due to 
temperature variations and thermal shock, and other effects such as chemical reaction or 
radiation attack. They might also include other time-dependent primitive variables such as 
creep, mechanical fatigue, thermal aging, thermal fatigue, or even effects such as seasonal 
attack (see Appendix A, Primitive Variables, Symbols, and Units). For most of these 
primitive variables, strength has been observed to decrease with an increase in the variable. 

The postulated constitutive equation accounts for the degradation of strength due to 
these primitive variables. The general form of the equation is 


! n 
-=n 

o i=l 


_VA 

A iF‘ A iO 


a. 

* 


( 1 ) 


where Aj, Aip and Ajo are the current, ultimate and reference values of a particular effect, 
ai is the value of an empirical constant for the i* effect or primitive variable, n is the 
number of product terms of primitive variables in the model, and S and So arc the current 
and reference values of material strength. Each term has the property that if the current 
value equals the ultimate value, the current strength will be zero. Also, if the current value 
equals the reference value, the term equals one and strength is not affected by that variable. 


This deterministic constitutive model may be calibrated by an appropriately curve-fitted 
least squares multiple linear regression of experimental data [7], perhaps supplemented by 
expert opinion. Ideally, experimental data giving the relationship between effects and 
strength is obtained. For example, data for just one effect could be plotted on log-log 
paper. A good fit for the data is then obtained by a linear regression analysis. This is 
illustrated schematically in Figure 1. The postulated constitutive equation, for a single 


effect, is then obtained by noting the linear relation between log S and 


log 


Af- Aq 


L Ap - A J 
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Fig. 1 Schematic of experimental data illustrating 

the effect of one primitive variable on strength. 



as follows: 


rvAoi 

log S - - 1 log I A-Xj+JosS 0 

iwi 

log S - log S Q > - a log I 




( 2 ) 


Note that the above equation (2) is for a primitive variable that lowers strength. If a 
variable raises strength, the exponent is negative. 

This general constitutive model may be used to estimate the strength of an aerospace 
propulsion system component under the influence of a number of diverse effects or 
primitive variables. The probabilistic treatment of this equation includes randomizing the 
deterministic multifactor interaction constitutive equation, performing probabilistic analysis 
by simulation and generating probability density function (p.d.f.) estimates for strength 
using a non-parametric method, maximum penalized likelihood [8,9]. Integration yields 
the cumulative distribution function (c.cLf.) from which probability statements regarding 
strength may be made. This probabilistic constitutive model predicts the random strength 
of an aerospace propulsion component due to a number of diverse random effects. 

This probabilistic constitutive model is embodied in two FORTRAN programs, 
PROMISC (Pro babilistic Material Strength Calibrator) and PROMISS (Prob abilistic 
Material Strength Simulator); see Final Technical Report, APPENDIX 1. PROMISS 
calculates the random strength of an aerospace propulsion component due to as many as 
eighteen diverse random effects. Results arc presented in the form of probability density 
functions and cumulative distribution functions of normalized strength, S/So- PROMISC 
calculates the values of the empirical material constants, a*. 

PROMISS (see Final Technical Repeat, APPENDIX 1) includes a relatively simple 
"fixed" model as well as a "flexible" model. The fixed model postulates a probabilistic 
constitutive equation that considers the primitive variables given in Table 1. The general 
form of this constitutive equation is given in equation (1), wherein there are now n = 7 
product terms, one for each effect or primitive variable listed above. Note that since this 
model has seven primitive variables, each containing four values of the variable, it has a 
total of twenty-eight variables. The flexible model postulates a probabilistic constitutive 
equation that considers up to as many as n = 18 product terms for primitive variables. 

These variables may be selected to utilize the theory and experimental data currently 
available for the specific strength degradation mechanisms of interest. The specific effects 
included in the flexible model are listed in Table 2. Note that in order to provide for future 
expansion and customization of the flexible model, six "other" effects have been provided. 
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Table 1 Primitive variables available in the fixed model 


i*ti Primitive 
^^ariable^^ 

Primitive 
Variable Type 

1 

Stress due to static load 

2 

Temperature 

3 

Chemical reaction 

4 

Stress due to impact 

5 

Mechanical fatigue 

6 

Thermal fatigue 

7 

Creep 

Table 2 Primitive variables available in the flexible model 


A. Environmental Effects 

1. Mechanical 

a. Stress 

b. Impact 

c. Other Mechanical Effect 

2. Thermal 

a. Temperature Variation 

b. Thermal Shock 

c. Other Thermal Effect 

3. Other Environmental Effects 

a. Chemical Reaction 

b. Radiation Attack 

c. Other Environmental Effect 

B. Time-Dependent Effects 

1. Mechanical 

a. Creep 

b. Mechanical Fatigue 

c. Other Mech. Time-Dep. Effect 

2. Thermal 

a. Thermal Aging 

b. Thermal Fatigue 

c. Other Thermal Time-Dep. Effect 

3. Other Time -Dependent Effects 

a. Corrosion 

b. Seasonal Attack 

c. Other Time-Dep. Effect 
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The considerable scatter of experimental data and the lack of an exact description of 
the underlying physical processes for the combined mechanisms of fatigue, creep, 
temperature variations and so on, make it natural, if not necessary to consider probabilistic 
models for a strength degradation model. Therefore, the fixed and flexible models 
corresponding to equation (1) are "randomized", and yield the "random normalized material 
strength due to a number of diverse random effects or primitive variables." Note that for 
the fixed model, equation (1) has the following form: 


S/So = f(AiF. Ai, Aio, A2F, A2, A20,-, A7F, A7, A70) (3) 


where Ai, Ajf and Ajo are the ultimate, current and reference values of the i* of seven 
effects or primitive variables as given in Table 1. In general, this expression can be written 
as, 


S/So = f(Xi), i = 1,..., 28, (4) 


where the Xj are the twenty eight independent variables in equation (3). Thus the fixed 
model is "randomized” by assuming all the independent variables, Xi, i = 1,..., 28, to be 
random and stochastically independent. For the flexible model, equation (1) has a form 
analogous to equations (3) and (4), except that there are as many as seventy-two 
independent variables. Applying probabilistic analysis to either of these randomized 
equations yields the distribution of the dependent random variable, normalized material 
strength, S/So- 

Although a number of methods of probabilistic analysis are available, [8] simulation 
was chosen for PROMISS. Simulation utilizes a theoretical sample generated by numerical 
techniques for each of the independent random variables. One value from each sample is 
substituted into the functional relationship, equation (3), and one realization of normalized 
strength, S/So, is calculated. This calculation is repeated for each value in the set of 
samples, yielding a distribution of different values for normalized strength. 

A probability distribution function is generated from these different values of 
normalized strength using a non-parametric method, maximum penalized likelihood. 
Maximum penalized likelihood generates the p.df. estimate using the method of maximum 
likelihood together with a penalty function to smooth it [9]. Finally, integration of the 
generated p.d.f. results in the cumulative distribution function, from which probabilities of 
normalized strength can be directly observed. 

PROMISS includes computational algorithms for both the fixed and the flexible 
probabilistic constitutive models. As described above, PROMISS randomizes the 
following equation: 


s n 

f-n 

*0 i=i 


.-A- 


A- C * A;. 


(5) 
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where 




a. 

. i 


Aj-Ajo 


is the i* effect, Ajf, Aio and Ao are random variables, a* is the i 0 * empirical material 
constant and S/So is normalized strength. There are a m a ximum of eighteen possible 
effects or primitive variables that may be included in the model For the flexible model 
option, they may be chosen by the user from those in Table 2. For the fixed model option, 
the primitive variables of Table 1 arc chosen. Within each primitive variable term the 
current, ultimate and reference values and the empirical material constant may be modeled 
as either deterministic (empirical, calculated by PROMISC), normal, lognormal or Wiebull 
random variables. Simulation is used to generate a set of realizations for normalized 
random strength, S/So, fro™ a set of realizations for primitive variables and empirical 
material constants. Maximum penalized likelihood is used to generate an estimate for the 
p.d.f. of normalized strength, from a set of realizations of normalized strength. Integration 
of the p.d.f. yields the c.d.f. Plot files are produced to plot both the p.d.f. and the c.d-f. 
PROMISS also provides information on S/SO statistics (mean, variance, standard deviation 
and coefficient of variation). A resident database, for database rather than user input of 
empirical material constants, is also provided. 

PROMISC performs a multiple linear regression on experimental data for as many 
as eighteen effects or primitive variables, yielding regression coefficients that are the 
empirical material constants, ai, required by PROMISS. 

It produces the multiple linear regression of the log transformation of equation (3), the 
PROMISS equation. When transformed it becomes 



( 6 ) 


or 


log S = log S 0 + £* log [-AjEcA.1 , 
m LAiF - AioJ 


(7) 


where 


A*' A V 


Aj:' Ajo 


is the i^ effect, Ai, Ajp and Ajo are primitive variable data and ai is the i^ 1 empirical 
material constant, or the i 1 * 1 regression coefficient to be predicted by PROMISC. Also, 
log So is the log transformed reference value of strength , or the intercept regression 
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coefficient to be predicted by PROMISC, and log S is the log transformed strength. 
Experimental data for up to eighteen possible effects, as given in Table 2, may be included. 
The primitive variable data may be either actual experimental data or expert opinion, directly 
read from input, or simulated data where expert opinion is specified as the mean and 
standard deviation of a normal or lognormal distribution. The simulated data option for 
input data was used in the early stages of code development to verify correct performance. 
The input data, whether actual or simulated, is read in and assembled into a data matrix. 
From this data matrix, a corrected sums of squares and crossproducts matrix is computed. 
From this sums of squares and crossproducts matrix, and a least squares methodology, a 
multiple linear regression is performed to calculate estimates for the empirical mater ial 
constant, a;, and the reference strength. So- These are the regression coefficients. 

PROMISC includes enhancements of the multiple linear regression analysis to 
screen data from "outliers" and coUinearities, determine "how well" the data fit the 
regression, quantify the importance and relative importance of each factor in die postulated 
constitutive equation, eq. (1), as well as check assumptions inherent in the use of multiple 
linear regression. Further details are provided in the Final Technical Report, Section 6.0, 
NASA Grant No. NAG 3-867, Supp. 2, "Probabilistic Lifetime Strength of Aerospace 
Materials via Computational Simulation." 
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3.0 SAMPLE PROBLEMS AND INPUT DATA 


Data input for PROMISC is user friendly and easy to enter. Actual experimental 
data can be input by the user or the user can select simulated experimental d ata . Two 
examples follow (see also Section 6.0, Appendix B). 


3.1 PROMISC Input For Fixed Model With Actual Experimental Data (40 Data Points) 

The 1st line of input (format 2E12.4, see item 1 below) determines the random 
number generator seed* and the data sample size. The 2nd line (format 13) determines 
either a fixed or a flexible model. The 3rd line (format 13) chooses the dependent variable 
as strength. The 4th line (format I3,2X,I3,2X,I3) determines the source of the data (actual 
experimental, normal or lognormal simulated data). The 5th through the 28th lines (format 
5D12.4) specify the data for the primitive variables. This sequence of twenty-five lines 
repeats for the other primitive variables. Lastly, the 179th through the 187th lines specify 
the data for the strength variable. A table listing the primitive variables, their units and 
symbols is given in Section 5.0, Appendix A. IMSL parameters are entered as indicated in 
items 12 and 13. 

1. Line 1 selects the Random Number Generator Seed (ISEED) and Sample Size (NTOT) 
EXAMPLE: 

12345 67 8 9012345678 90** 

_ 4 o 


2. Line 2 selects either Fixed or Flexible Model (MODEL = 0 is flexible, MODEL - 1 is 
fixed). 

EXAMPLE: 

12345678901234567890 

1 

3. Line 3 selects strength as the dependent variable (DEPV = 1 is strength dependent 
variable, DEPV = 0 is not strength dependent variable). 

EXAMPLE: 

12345678901234567890 

1 


* If actual experimental data are used (AFDS = 0; ADS — 0; AODS = 0) the random number generator seed, 
ISEED, is read in but not used in the program. 

** NOTE: the ruler is to aid the user in formatting and is not a part of the input. 
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4a. Line 4 specifies the source of the data for the three variables within the Quasi-static 
Stress Effect by using flags. The flag names are AFDS, ADS and AODS. Setting a flag to 
0 indicates actual experimental data is directly read from input Setting a flag to 1 indicates 
simulated experimental data, normally distributed, is generated by the program via a 
random number generator. Setting a flag to 2 indicates simulated experimental data, 
lognormally distributed, is generated by the program via a random number generator. 
Setting a flag to 3 indicates data not available. 

EXAMPLE: 

12345678901234567890 
0 0 0 


4b. Lines 5 to 28 contain the actual experimental data for the three variables within the 
Quasi-static Stress Effect The Quasi-static Stress Effect variable names have FORTRAN 
names given in Section 5.0, APPENDIX A. 

EXAMPLE: 


123456789012345678901234567890123456789012345678901234567890 


0 . 1179D+03 
0 . 1292D+03 
0 . 1321D+03 
0 . 1278D+03 


0 . 1344D+03 0 . 1273D+03 0.1260D+03 
0 . 1295D+03 0 . 1241D+03 0.1265D+03 
0 . 1387D+03 0 . 1370D+03 031264D+03 
0 . 1269D+03 etc. 


0 . 1403D+03 
0.1229D+03 
0. 1323D+03 


(All of the values for this and other effects are too numerous to include in this description. 
See Section 6.0 APPENDIX B.) 


5. Lines 29 to 53 input data for the Impact Effect in the same manner as was described for 
the Quasi-static Stress Effect The Impact Effect variable names have FORTRAN names 
given in Section 5.0, APPENDIX A. 

EXAMPLE: 


12345678901234567890 


0 0 0 
0 . 9870D+00 
0.9977D+00 
0 . 9991D+00 
0 . 1000D+00 


0 . 9966D+00 
0 . 9950D+00 
0 . 1000D+01 
0 . 1001D+01 


0 . 1002D+01 
0 . 1001D+01 
0 . 1003D+01 


0 . 9997D+00 
0 . 1001D+01 
0 . 9945D+00 


etc. 


0 . 1000D+01 
0 . 1005D+01 
0 . 9952D+00 
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6. Lines 54 to 78 input data for the Temperature Effect in the same manner as was 
described for the Quasi-static Stress Effect. The Temperature Effect variable names have 
FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 


i 9^am;7RQ0i 7 345 67 8 901 ? 34 5 67 8 9012 34 5 67 8 901234 5 678901 2 34 5 678 20, 


0 0 0 
0 .2378D+04 
0 . 2668D+04 
0 . 2708D+04 
0 .2738D+04 


0 . 2640D+04 
0 . 2598D+04 
0 .2739D+04 
0.2768D+04 


0 . 2789D+04 
0 .2770D+04 
0 . 2811D+04 
etc. 


0.2723D+04 
0 . 2770D+04 
0 . 2583D+04 


0 . 2739D+04 
0.2856D+04 
0.2600D+04 


7 Lines 79 to 103 input data for the Chemical Effect in the same manner as was described 
for the Quasi-static Stress Effect. The Chemical Effect variable names have FORTRAN 
names given in Section 5.0, APPENDIX A. 

EXAMPLE: 


1 7 R 4 . 5 67 8 901234 5 67 8 90 1234 5 678 901234 5678 9Q1234 5 678 90 1234 5 678 9Q 


0 0 0 
0 . 9870D+00 
0 . 9977D+00 
0 .9991D+00 
0 . 1000D + 01 


0 . 9966D+00 
0 . 9950D+00 
0 . 1000D+01 
0 . 1001D+01 


0 . 1002D+00 
0 . 1001D+01 
0 . 1003D+01 
etc. 


0 . 9997D+00 
0 . 1001D+01 
0 . 9945D+00 


0 . 1000D+01 
0 . 1005D+01 
0.9952D+00 


8 Lines 104 to 128 input data for the Creep Effect in the same manner as was described for 
the Quasi-static Stress Effect. The Creep Effect variable names have FORTRAN names 
given in Section 5.0, APPENDIX A. 


EXAMPLE: 


19 34 5678 9012345 67890123456 7 8901234 5678901234567890 12 34 5 678 90 


0 0 0 
0 . 9072D+04 
0 . 9936D+04 
0 . 101 6D+05 
0 . 9827D+04 


0 .1034D+05 
0 . 9959D+04 
0 . 1067D+05 
0 . 9758D+04 


0 . 9792D+04 
0. 9544D+04 
0 . 1054D+05 
etc. 


0 . 9696D+04 
0.9733D+04 
0 . 9723D+04 


0 . 1080D+05 
0 . 9450D+04 
0 . 1018D+05 
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9. Lines 129 to 153 input data for the Impact Effect in the same manner as was described 
for the Quasi-static Stress Effect The Impact Effect variable names have FORTRAN 
names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

12345 67 8 90 12 345 678 90 12345678 90123456789 0 12345 67 890123 45 67890 
0 0 0 

0 . 5749D+01 0.7465D+01 0.6696D+01 0.6565D+01 0.8136D+01 

0 . 6894D+01 0 . 6926D+01 0.6362D+01 0.6615D+01 0.6237D+01 

0 . 7214D+01 0 . 7947D+01 0.7752D+01 0.6602D+01 0.7235D+01 

0 . 6744D+01 0. 6649D+01 etc. 


10. Lines 154 to 178 input data for the Thermal Fatigue Effect in the same manner as was 
described for the Quasi-static Stress Effect The Thermal Fatigue Effect variable names 
have FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 


12345678901234567890 1234567890123456789012345678901 234567890 
0 0 0 

0 . 2464D+01 0.3199D+01 0.2870D+01 
0 . 2954D+01 0.2968D+01 0.2727D+01 
0 . 3092D+01 0 . 3406D+01 0.3322D+01 
0 . 2890D+01 0 . 2850D+01 etc. 


0.2813D+01 
0 . 2835D+01 
0 . 2830D+01 


0 . 3487D+01 
0 . 2673D+01 
0 . 3101D+01 


11a. Line 179 specifies the source of the strength data by using flags. The flag names are 
AFDS, ADS and AODS. Setting a flag to 0 indicates actual experimental data is directly 
read from input Setting a flag to 1 indicates simulated experimental data, normally 
distributed, is generated by the program via a random number generator. Setting a flag to 2 
indicates simulated experimental data, lognormally distributed, is generated by the program 
via a random number generator. Setting a flag to 3 indicates data not available. 

EXAMPLE: 


123456789Q 

0 

1 lb. Lines 180 to 187 contain the actual experimental data for strength. The FORTRAN 
name for strength is STR. 

EXAMPLE: 


123456789012345678901234567890123456789012345678901234567890 


0 . 1117D+02 
0 .0831D+02 
0 . 0798D+02 
0 . 0774D+02 


0 . 0855D+02 
0 .0895D+02 
0 . 0776D+02 
0 . 0751D+02 


0 . 0735D+02 
0 . 0749D+02 
0 . 0719D+02 
etc. 


0 . 0786D+02 
0 . 074 9D+02 
0 . 0906D+02 


0 . 0774D+02 
0 . 0687D+02 
0 . 0891D+02 
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12. The CORVC [1] parameters are IDO, LDX, EFRQ, IWT, MOPT, I COPT, LDCOV, 
and LDINCD and are entered in that order as follows (line 188): 

EXAMPLE: 

I934 5fi7fl9m 234 5 678 9 012345678901234567891) 

0 240 0 0 0 1 19 19 


13. The RCOV [1] parameters are INTCEP, NOND, NDEP, LDCOV, LDB, LSR, and 
LDSCPE and are entered in that order as follows (line 189): 

EXAMPLE: 

1 2745 67890123456789012345678 9012345 67 890 
1 7 1 19 19 19 1 


3.1.1 Discussion of Results 

Execution of PROMISC (source code entitled PROMISC5.FOR) produces an 
output file that gives computed values for the empirical material constants, aj. Output also 
includes statistics and diagnostics that enhance the PROMISC multiple linear regression 
analysis (see Section 2.0, Theoretical Background and/or Final Technical Report, Section 6 
Increased Capability of PROMISC, NASA Grant No. NAG 3-867, Supp. 2, "Probabilistic 
Lifetime Strength of Aerospace Materials via Computational Simulation"). 

Execution of the PROMISC code for the example problem whose input is specified 
in Section 3.1 produces an output file that includes regression or least squares estimates of 
the exponents (empirical material constants, ai) in the multifactor interaction equation under 
the heading BBBB. For example, the constant coefficient (y-intcrcept, log So) is predicted 
to be 0.6717, the exponent for quasi-static stress to be 0.6449, the exponent for impact to 
be 0.0000, and so on. 

To determine "how well" the data fit the regression equation consider the next 
subheading in the output file, where R-squared is given. R-squared and the Adjusted R- 
squared columns give the percentage of the variation in the data which can be explained by 
the multiple linear regression equation. In this example the multiple linear regression 
equation explains about 99.7% of the variation. An R-squared value of 100% indicates an 
exact fit, whereas an R-squared value of 0 % indicates no correlation between the 
independent variable (strength) and any of the regression variables (effects or primitive 
variables). This is a first assessment of the adequacy of the model. 

In the next heading under ANOVA Table, an alternative statistical measure of "how 
well" the data fit the "complete" multiple linear regression equation is given. This statistical 
measure is given in the last column headed by "Prob. of Larger F". This term, known as a 
p- value, gives the probability that the relationship was due to chance alone and not to the 
multiple linear regression equation. In this case the probability is virtually zero, which 
corroborates the observations from the R-squared columns. P-values near zero indicate 
that the multifactor interaction equation is significant Thus, we conclude that the multiple 
linear regression equation is quite effective in predicting strength. 
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Among the seven regression variables (effects or primitive variables) it is desirable 
to assess the relative importance of each. From the section headed by BBBB or the one 
headed by "Inference on Coefficients", the estimates of the exponents in the multifactor 
interaction equation can be viewed. Since the exponent for impact is virtually zero it can be 
reasonably assumed that this factor does little to predict the strength. However, it is not 
wise to be mislead by the mere size of the coefficients. The statistical terms which one 
should focus on are the "F-statistic" column and the "Prob. of Larger F" column in the 
section headed by Sequential Statistics. For the first variable, quasi-static stress, the p- 
value gives the probability that the factor is not significant in the prediction equation. The 
exponent for quasi-static stress is 0.6449 and the associated p-value is 0.0660, which 
indicates quasi-static stress is significant In most statistical circles, p-values near 0.05 or 
less are considered significant Moreover, the sixth effect mechanical fatigue, has a 
coefficient (exponent of -0. 1248) but its p-value is virtually zero. This result indicates that 
mechanical fatigue (because of its smaller p-value), has greater predictive capability than the 
quasi-static stress effect However, both should be considered significant The remaining 
factors are completely insignificant. 

In outlier detection, use the "Case Analysis" section and the factor, "Cook's D". 
Cook’s D is given first in the second line of each of the forty cases. For example, for 
observation 1 Cook's D is 37.7568. The larger the value of Cook's D, the farther the 
observation is away from the center of the data. We consequently note that observation 14, 
with a Cook's D of 20.0449, has the smallest value, and observation 34, with Cook’s D of 
50.3682, has the largest value. However even the largest value is well within acceptable 
ranges indicating that there are no outliers in the data. 

The residual column and the standardized residual column can be used to test the 
hypothesis that errors are normally distributed. In addition, the remaining columns can be 
used for many statistical purposes in verifying the assumptions in multiple linear 
regression. 


92 



3.2 PROMISC Input For Fixed Model With Actual Experimental Data (240 Data Points) 

The 1st line of input (format 2E12.4, see item 1 below) determines the random 
number generator seed* and the data sample size. The 2nd line (format 13) determines 
either a fixed or a flexible model. The 3rd line (format 13) chooses the dependent variable 
as strength. The 4th line (format I3,2X,I3,2X,I3) determines the source of the data (actual 
experimental, normal or lognormal simulated data). The 5th through the 148m lines 
(format 5D12.4) specify the data for the primitive variables. This sequence of one hundred 
forty-five lines repeats for the other primitive variables. Lastly, the 1019th through the 
1067th lines specify the data for the strength variable. A table listing the primitive 
variables, their units and symbols is given in Section 5.0, Appendix A. IMSL parameters 
are entered as indicated in items 12 and 13. 

1. Line 1 selects the Random Number Generator Seed (ISEED) and Sample Size (NTOT). 
EXAMPLE: 

1234 5 67 8901234567890** 

1 240 — — ~ 


2. Line 2 selects either Fixed or Flexible Model (MODEL = 0 is flexible, MODEL = 1 is 
fixed). 

EXAMPLE: 

12345678901234567890 

1 


3. Line 3 selects strength as the dependent variable (DEPV — 1 is strength dependent 
variable, DEPV = 0 is not strength dependent variable). 

EXAMPLE: 

12345678901234567890 

1 


* If ac tual experimental data are used (AFDS = 0; ADS = 0; AODS = 0) the random number generator seed, 
ISEED, is read in but not used in the program. 

** NOTE: the ruler is to aid the user in formatting and is not a part of the input. 


93 


4a. Line 4 specifies the source of the data for die three variables within the Quasi-stadc 
Stress Effect by using flags. The flag names are AFDS, ADS and AODS. Setting a flag to 
0 indicates actual experimental data is directly read from input Setting a flag to 1 indicates 
simulated experimental data, normally distributed, is generated by the program via a 
random number generator. Setting a flag to 2 indicates simulated experimental data, 
lognonnally distributed, is generated by the program via a random number generator. 
Setting a flag to 3 indicates data not available. 

EXAMPLE: 

12345678901234567890 
0 0 0 


4b. Lines 5 to 148 contain the actual experimental data for the three variables within the 
Quasi-static Stress Effect. The Quasi-static Stress Effect variable names have FORTRAN 
names given in Section 5.0, APPENDIX A. 

EXAMPLE: 


123456789012345678901234567890123456789012345678901234567890 


0.1179D+03 
0 . 1292D+03 
0 . 1321D+03 
0 . 1278D+03 


0.1344D+03 0 . 1273D+03 0.1260D+03 
0 . 1295D+03 0 . 1241D+03 0.1265D+03 
0 . 1387D+03 0 . 1370D+03 031264D+03 
0.1269D+03 etc. 


0 . 1403D+03 
0 . 1229D+03 
0 . 1323D+03 


(All of the values for this and other effects are too numerous to include in this description. 
See Section 6.0 APPENDIX B.) 


5. Lines 149 to 293 input data for the Impact Effect in the same manner as was described 
for the Quasi-static Stress Effect The Impact Effect variable names have FORTRAN 
names given in Section 5.0, APPENDIX A. 

EXAMPLE: 


123 4 5678901? 34567890 


0 0 0 
0 . 9870D+00 
0 . 9977D+00 
0 . 9991D+00 
0 . 1000D+00 


0 . 9966D+00 
0 . 9950D+00 
0 . 1000D+01 
0 . 1001D+01 


0 . 1002D+01 
0 . 1001D+01 
0 . 1003D+01 


0 . 9997D+00 
0 . 1001D+01 
0 . 9945D+00 


etc. 


0 . 1000D+01 
0 . 1005D+01 
0. 9952D+00 
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0 0 0 

0.9072D+04 0 . 1034D+05 0.9792D+04 0.9696D+04 0.1080D+05 
0 . 9936D+04 0.9959D+04 0.9544D+04 0.9733D+04 0.9450D+04 
0 . 1016D+05 0 . 1067D+05 0.1054D+05 0.9723D+04 0.1018D+05 
0 . 9827D+04 0.9758D+04 CtC. 


95 


9. Lines 729 to 873 input data for the Impact Effect in the same manner as was described 
for the Quasi-static Stress Effect The Impact Effect variable names have FORTRAN 
names given in Section 5.0, APPENDIX A. 

EXAMPLE: 

1234 5 67 89012345 67 8 90 12345678 90 1234 5 67 8 90 12345 67 8 90 12345 67 890 
0 0 0 

0 . 5749D+01 0.7465D+01 0.6696D+01 0.6565D+01 0.8136D+01 

0 . 6894D+01 0 . 6926D+01 0.6362D+01 0.6615D+01 0.6237D+01 

0 . 7214D+01 0 . 7947D+01 0.7752D+01 0.6602D+01 0.7235D+01 

0 . 6744D+01 0 . 6649D+01 etc. 


10. Lines 874 to 101 8 input data for the Thermal Fatigue Effect in the same manner as was 
described for the Quasi-static Stress Effect The Thermal Fatigue Effect variable names 
have FORTRAN names given in Section 5.0, APPENDIX A. 

EXAMPLE: 


123456789012345678901234567890123456789012345678901234567890 


0 

0, 

0, 

0 . 

0 , 


0 0 
24 64D+01 
2954D+01 
3092D+01 
2890D+01 


0 . 3199D+01 0.2870D+01 
0.2968D+01 0 . 2727D+01 
0 . 3406D+01 0 . 3322D+01 
0 . 2850D+01 etc. 


0 .28130+01 0 
0 . 2835D+01 0, 
0 . 2830D+01 0 


3487D+01 

2673D+01 

3101D+01 


1 la. Line 1019 specifies the source of the strength data by using flags. The flag names are 
AFDS, ADS and AODS. Setting a flag to 0 indicates actual experimental data is directly 
read from input Setting a flag to 1 indicates simulated experimental data, normally 
distributed, is generated by the program via a random number generator. Setting a flag to 2 
indicates simulated experimental data, lognormally distributed, is generated by the program 
via a random number generator. Setting a flag to 3 indicates data not available. 

EXAMPLE: 


1234567890 

0 

1 lb. Lines 1020 to 1067 contain the actual experimental data for strength. The 
FORTRAN name for strength is STR. 

EXAMPLE: 


12345 678 9Q12345 678 9Q123456789Q12345 678 9Q12345 678 9Q12345 678 90 


0 .1112D+02 
0 . 8277D+01 
0 . 7954D+01 
0 .7711D+01 


0 .8518D+01 
0 .8921D+01 
0 . 7729D+01 
0 . 7478D+01 


0 . 7325D+01 
0 . 7464D+01 
0 . 7166D+01 
etc. 


0.7826D+01 
0 . 7464D+01 
0 . 9028D+01 


0 . 7708D+01 
0 . 6842D+01 
0.8875D+01 
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12. The CORVC [1] parameters are IDO, LDX, IFRQ, IWT, MOPT, ICOPT, LDCOV, 
and LDINCD and are entered in that order as follows (line 1068): 

EXAMPLE: 

1 ?34567BQCn 2345 6 78901234 56789012345678 90 
0 240 0 0 0 1 19 19 


13. The RCOV [1] parameters are INTCEP, NOND, NDEP, LDCOV, LDB, LSR, and 
LDSCPE and are entered in that order as follows (line 1069): 

EXAMPLE: 

1734567890123456789012345678 901234567890 
1 7 1 19 19 19 1 


3.2. 1 Discussion of Results 

Execution of PROMISC (source code entitled PROMISC5.FOR) produces an 
output file that gives computed values for the empirical material constants, ai. Output also 
includes statistics and diagnostics that enhance the PROMISC multiple linear regression 
analysis (see Section 2.0, Theoretical Background and/or Final Technical Report, Section 6 
Increased Capability of PROMISC, NASA Grant No. NAG 3-867, Supp. 2, "Probabilistic 
Lifetime Strength of Aerospace Materials via Computational Simulation"). 

Execution of the PROMISC code for the example problem whose input is specified 
in Section 3.2 produces an output file that includes regression or least squares estimates of 
the exponents (empirical material constants, aj) in the multifactor interaction equation under 
the h eadin g BBBB. For example, the least squares estimate of the constant term (y- 
intercept, log So) is predicted to be 4.99, quasi-static stress to be 0. 14, mechanical fatigue 
to be -0.05, and so on. 

To see "how well" the data fit the regression equation, consider the subheading 
which contains R-squared. From the R-squared coefficients, it can be observed that about 
99.99% of the variation in the data can be explained by the regression equation. Recall that 
100% is the value when the model is exact. Under the next subheading, the ANOVA Table 
helps to quantify the model adequacy in terms of probability. Under the source 
"regression" it is observed that the p-value is virtually zero. The asterisks in the 
"Overall F" column indicate the F-statistic exceeded 99,999.99 and therefore could not be 
displayed in the space allocated. 

Among the seven factors (effects or primitive variables) in the multifactor 
interaction equation, we can associate relative importance with the size of the p-values in 
the "Sequential Statistics" Table. The smaller the p-value (i.c., "Prob. of Larger F ) then 
the more significance associated with the factor. Notice that factors 1, 2 and 5 are 
significant with p-values below the 0.05 level. However factors 3, 4, 6 and 7 are 
insignificant. Therefore, it may be beneficial to reanalyze the data by deleting factors 3, 4, 
6 and 7. 
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To look for outlying observations, the "Case Analysis" Table is used. The column 
headed by "Cook's D” is used for the detection of aberrant observations. The output is 
complicated in this case by the number of observations present in the data set As 
explained in Section 3.1.1, the larger values of "Cook's D" indicate observations which are 
further away from the center of the data. However, when 240 observations are present the 
process of scanning the values of Cook's D is more tedious. Observation 85 is the one 
with the largest Cook's D of 0.0592 and standardized residual of 2.6871, both of which 
are within acceptable limits. 

The residuals and standardized residuals can be used to test the assumptions of 
normality for the errors. 
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3.3 General PROMISC Program Notes 


1 . 

2 . 


Any effect data, \t- not included in the model is set to zero. 

Fatigue cycles should be input as log cycles rather than as cycles. Reasons are the 
same as for PROMISS. 


3. Strength response, S, should be input as strength in MPa or ksi. PROMISC 

will internally calculate log S. 

4 . The stress data effect in PROMISC has not been customized for negative 
lognormal values as it has been for PROMISS. 
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5.0 APPENDIX A 

PRIMITIVE VARIABLES, SYMBOLS, AND UNITS 
Table A 1.2 Primitive variables, symbols, and units for PROMISS 


Primitive Variables 
(Effect) 

Theory 

Symbol 

FORTRAN 

Name 

Units 

SI 

U.S 

OUASI-STATIC 
STRESS EFFECT 

Ultimate value 

Ssf 

AF 

MPa 

ksi 

Current value 

OS 

A 

MPa 

ksi 

Reference value 

oso 

AO 

MPa 

ksi 

IMPACT EFFECT 

Ultimate value 

Sdf 

AF2 

Dimensionless 


Current value 

gd 

A2 

Dimensionless 


Reference value 

ODO 

A02 

Dimensionless 


OTHER 

MECHANICAL 

EFFECTS 

Ultimate value 

A3F 

AF3 

Dimensionless 


Current value 

A 3 

A3 

Dimensionless 


Reference value 

A 30 

A03 

Dimensionless 


TEMPERATURE 

Ultimate value 

t f 

AF4 

°C 

op 

Current value 

T 

A4 

°C 

OF 

Reference value 

T 0 

A04 

°C 

°F 

THERMAL 

SHOCK 

Ultimate value 

Asp 

AF5 

Dimensionless 


Current value 

As 

A5 

Dimensionless 


Reference value 

A 50 

A05 

Dimensionless 
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Primitive Variables Theory 
(Effect) Symbol 


FORTRAN 

Name 


SI 


Units 


U.S. 


OTHER 

THERMAL 

EFEEC1S. 

Ultimate value 

A6F 

AF6 

Dimensionless 


Current value 

A6 

A6 

Dimensionless 


Reference value 

A<50 

A06 

Dimensionless 


CHEMICAL 

REACTION 

Ultimate value 

Rf 

AF7 

Dimensionless 


Current value 

R 

A7 

Dimensionless 


Reference value 

Ro 

A07 

Dimensionless 


RADIATION 

Ultimate value 

AgF 

AF8 

Dimensionless 


Current value 

As 

A8 

Dimensionless 


Reference value 

Ago 

AO 8 

Dimensionless 


OTHER EFFECTS 

Ultimate value 

A 9 F 

AF9 

Dimensionless 


Current value 

A 9 

A9 

Dimensionless 


Reference value 

A 90 

AQ9 

Dimensionless 


CREEP 

Ultimate value 

tCF 

AF10 

Hours 

Hours 

Current value 

tc 

A10 

Hours 

Hours 

Reference value 

kx) 

AO10 

Hours 

Hours 

MECHANICAL 

FATIGUE 

Ultimate value 

Nmf 

AF11 

log of cycles 


Current value 

Nm 

All 

log of cycles 


Reference value 

NmO 

AOll 

log of cycles 
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Primitive Variables Theory 
(Effect) Symbol 


FORTRAN 

Name 


SI 


Units 


U.S. 


OTHER TIME 
DEPENDENT 
MECHANICAL 
EFFECT 



Ultimate value 

Ai2F 

AF12 

Current value 

A 12 

A12 

Reference value 

A 120 

A012 

THERMAL AGING 



Ultimate value 

Ai3F 

AF13 

Current value 

A 13 

A13 

Reference value 

Ai30 

A013 

THERMAL 

FATIGUE 



Ultimate value 

Ntf 

AF14 

Current value 

Ny 

A14 

Reference value 

Nyo 

A014 

OTHER TIME- 
DEPENDENT 
THERMAL 
EFFECTS 



Ultimate value 

Al5F 

AF15 

Current value 

Al5 

A15 

Reference value 

Al50 

A015 

CORROSION 



Ultimate value 

Ai6F 

AF16 

Current value 

Ai6 

A16 

Reference value 

A 160 

A016 

SEASONAL 

ATTACK 



Ultimate value 

Ai7F 

AF17 

Current value 

A 17 

A17 

Reference value 

Ai70 

A017 
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Dimensionless 

Dimensionless 

Dimensionless 


Dimensionless 

Dimensionless 

Dimensionless 


log of cycles 
log of cycles 
log of cycles 


Dimensionless 

Dimensionless 

Dimensionless 


Dimensionless 

Dimensionless 

Dimensionless 


Dimensionless 

Dimensionless 

Dimensionless 


Primitive Variables 
(Effect) 


OTHER TIME- 
DEPENDENT 




Ultimate value 
Current value 
Reference value 


Theory FORTRAN Units 

Symbol Name SI U.S. 


Ai8F 

Ai8 

Ai80 


AF18 Dimensionless 

A18 Dimensionless 

A018 Dimensionless 
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6.0 APPENDIX B 

PROMISC SAMPLE PROBLEM: INPUT AND OUTPUT FILES 


Sample problems arc discussed in sections 3.1 and 3.2. One sample problem 
corresponds to each section. The input and output file names for each sample problem are 
listed below. The enclosed disk also includes die same files. 

3. 1 PROMISC Input For Fixed Model With Actual Experimental Data (40 Data Points) 

Input File(s): 31PRC4.INP 
Output File: 3 1 PRC4.0UT 

3.2 PROMISC Input For Fixed Model With Actual Experimental Data (240 Data Points) 

Input File(s): 32PRC5.INP 
Output File: 32PRC5.0UT 
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7.0 APPENDIX C 


1. RNSET - 

2. RNNOR - 

3. RNLNL - 

4. CORVC - 

5. RCOV - 

6. RSTAT - 

7. RCASE - 


IMSL SUBROUTINE CALLS FROM PROMISC 


Initializes a random seed for use in the IMSL 
random number generators. 

Generates pseudorandom numbers from a standard 
normal distribution using an inverse CDF 
method. 

Generates pseudorandom numbers from a 
lognormal distribution. 

Generates sums of squares and crossproducts. 

Performs the multiple linear regression. 

Computes statistics related to the regression. 

Computes further statistics and diagnostics related to the regression. 


8. Other IMSL Subroutines - SSCAL, SADD, UMACH, WRRRN, and WRIRN. 
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Figure A.1 COMPARISON OF VARIOUS LEVELS OF UNCERTAINTY OF 
STRESS ON (ksi) PROBABLE STRENGTH FOR A 
NICKEL-BASED SUPERALLOY AT ROOM TEMPERATURE 
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Figure A.2 COMPARISON OF VARIOUS LEVELS OF UNCERTAINTY OF 

MECHANICAL FATIGUE (CYCLES) ON PROBABLE STRENGTH 
FOR A NICKEL-BASED SUPERALLOY AT ROOM TEMPERATURE 
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Figure A.3 COMPARISON OF VARIOUS L£VELS OF UNCERTAINTY OF 

THERMAL FATIGUE (CYCLES) ON PROBABLE STRENGTH FOR 
A NICKEL-BASED SUPERALLOY AT ROOM TEMPERATURE 



figure A.4 COMPARISON OF VARIOUS LEVELS OF UNCEI 
OF CREEP (HOURS) ON PROBABLE STRENG1 
NICKEL-BASED SUPERALLOY AT ROOM TEMI 






Figure A.5 COMPARISON OF VARIOUS LEVELS OF UNCERTAINTY 
OF STRESS (ksi) ON PROBABLE STRENGTH FOR 
A NICKEL-BASED SUPERALLOY AT 0.57 Tm 



Figure A.6 COMPARISON OF VARIOUS LEVELS OF UNCERTAINTY OF 

MECHANICAL FATIGUE (CYCLES) ON PROBABLE STRENGTH 
FOR A NICKEL-BASED SUPERALLOY AT 0.57 Tm 
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Figure A.7 COMPARISON OF VARIOUS LEVELS OF UNCERTAINTY 
OF THERMAL FATIGUE (CYCLES) ON PROBABLE STRE 
FOR A NICKEL— BASED SUPERALLOY AT Tm 
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Figure A.8 COMPARISON OF VARIOUS LEVELS OF UNCERTAINTY 
OF CREEP (HOURS) ON PROBABLE STRENGTH 
FOR A NICKEL-BASED SUPERALLOY AT 0.57 Tm 


13.0 APPENDIX 4 : MATERIAL DATA PLOTS 


13.1 Introduction 

This appendix displays plots of material properties data collected from the open 
literature for certain aerospace materials. Included are tensile (quasi-static stress), fatigue 
and creep plots for INCONEL 718, tensile and fatigue plots for other mckel-base 
superalloys and creep data for metal matrix composites and reinforcement fibers. 


13.2 Tensile, Fatigue and Creep Data Plots of INCONEL 718 
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EFFECT OF TEMPERATURE ON TENSILE STRENGTH 

FOR INCONEL 71 8 
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EFFECT OF TEMPERATURE ON YIELD STRENGTH 
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EFFECT OF TEMPERATURE ON YIELD STRENGTH 

FOR INCONEL 718 



(!s>t) Hiowaais 


128 


TEMPERATURE (°F) 




TH IV H10N3H1S / (! s> l) H10N3U1S 

129 





EFFECT OF TEMPERATURE ON YIELD STRENGTH 
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EFFECT OF TEMPERATURE ON YIELD STRENGTH 

FOR INCONEL 71 8 
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EFFECT OF TEMPERATURE ON YIELD STRENGTH 

FOR INCONEL 71 8 
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Log (FCmax— FCmin)/(FCmax— FC) 



EFFECT OF CYCLES ON FATIGUE STRENGTH 
FOR INCONEL 718 
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EFFECT OF CYCLES ON FATIGUE STRENGTH 
FOR INCONEL 718 
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EFFECT OF CYCLES ON FATIGUE STRENGTH 
FOR INCONEL 718 
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Fatigue Cycles/Fat. Cyc. at Min Strength 
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EFFECT OF CYCLES ON FATIGUE STRENGTH 
FOR INCONEL 71 8 


li- 



(!S>|) H10N3dlS 3n0llVJ 


143 


FATIGUE CYCLES X 1000 



EFFECT OF CYCLES ON FATIGUE STRENGTH 
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Fatigue Cycles/Fat. Cyc. at Min Strength 


EFFECT OF CYCLES ON FATIGUE STRENGTH 
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EFFECT OF TIME ON RUPTURE STRENGTH 
FOR INCONEL 71 8 
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EFFECT OF TIME ON RUPTURE STRENGTH 
FOR INCONEL 718 


U- 



J 8> ! 'H1SN3&LS sarudna 


EFFECT OF TIME ON RUPTURE STRENGTH 
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